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ABSTRACT 


Peak  pressure  calibration,  pressure- time  profile  response  tests, 
and  pressure-resistance  hysteresis  of  the  Manganin  in  C- 7  epoxy 
transducer  are  described.  It  is  shown  that  predicted  and  measured 
pressure- time  profiles  agree  to  within  10  percent.  Also- presented 
are:  the  transformation  of  profile  measurements  from  transducer  to 

test  specimen,  Hugoniot  temperature  calculations  for  C-7  epoxy  up 
to  ~500  kbar  peak  pressure,  a  method  of  constructing  gages  of 
Manganin  in  conducting  materials,  and  gage  field  test  results  on 
the  FLAT  TOP  series  of  events. 
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I  INTRODUCTION 


A.  GENERAL 

Current  activity  in  underground  nuclear  detonations  has  generated 
considerable  interest  and  consequent  research  in  the  shock  propagation 
properties  of  rocks  and  soils.  Laboratory  measurements  typically  deter¬ 
mine  the  Hugoniot  characteristics  of  proposed  shot  media,  which  are  used 
in  field  test  predictions  of  stress-time  and  stress-distance  profiles. 
Field  instrumentation  has  been  developed  to  measure  the  latter,1  but  the 
former  has  been  difficult  to  achieve.  tn  fact,  stress-time  profile  mea¬ 
surements  of  explosively  produced  shock  waves,  even  on  a  laboratory  scale, 
have  been  limited  to  short  duration  and,  depending  on  the  type  of  trans¬ 
ducer,  to  rather  narrow  stress  ranges.2*3 

The  objective  of  the  investigation  discussed  in  this  report  is  the 
development  of  a  piezoresistive  transducer  or  system  of  transducers 
capable  of  measuring  stress-time  profiles  in  rocks  and  soils.  Specific 
response  goals  are  short  rise  time,  <0..!  fj.se c ,  long  duration  recording, 

>  50  /xsec,  and  an  overall  stress  accuracy  of  ±10  percent  in  a  stress 
range  from  10  kbar  to  as  high  as  possible,  preferably  at  least  500  kbar. 

Since  the  primary  goal  of  this  work  is  a  transducer  system  for  field 
measurements,  fundamental  investigations  of  the  piezoresistive  effect 
in  conductors  under  dynamic  compression  have  been  exa:  ined  only  with 
respect  to  the  goals  outlined  above.  Although  preliminary  studies  of 
resistivity  vs.  peak  pressure  have  been  made  for  several  alloys,  peak 
pressure  calibration,  pressure-resistance  hysteresis,  shock  heating,  etc. 
have  been  studied  for  only  one  type  of  alloy,  Dri ver-Harris  Manganin, * 
and  one  specific  wire  size,  0.003  inch  diameter.  Except  for  peak  pres¬ 
sure  calibration,  the  measurements  have  been  made  in  a  somewhat  cursory 
manner  in  that  only  a  few  experimental  tests  of  each  type  have  been  per¬ 
formed,  whereas  more  quantitative  results  would  require  a  statistical 
determination  and  much  more  elaborate  techniques.  These  would  result 
in  major  investigations  in  themselves. 
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It  should  be  noted  that  the  investigations  reported  Below  have  been 
performed  in  isotropic,  homogeneous  media  with  plane  shock  waves  aiid  « 

analyzed  using  one-dimensional  shock  wave  theory.  The  departure  from 
these  conditions  in  actual  field  tests  may  be  appreciable  due  to  geologi-  - 
cal  faults  and  inhomogeneities  in  the  region  of  transducer  placement. 

The  effect  of  these  unideal  co-ditions  on  transducer  response  can  be 
estimated,  at  present,  in  only  a  very  qualitative  manner  and  has  not 
been  undertaken  in  the  present  work.  Therefore,  field  test  results  in 
which  flow  in  the  vicinity  of  the  transducer  departs  greatly  from  one 
dimensionality  must  be  regarded  as  qualitative  also.  This  does  not 
include  spherical  divergence  from  a  “point  source”  detonation  since,  in 
general,  the  measurement  is  made  at  a  sufficiently  large  radius  so  that 
divergence  can  be  neglected. 

B.  BACKGROUND 

Prior  work  in  the  field  of  high  dynamic  stress  transducers  has  been 
reviewed  briefly  in  a  previous  report*  and  in  Reference  1.  Similarly, 
requirements  for  and  selection  of  piezoresistance  transducer  components 
have  been  outlined  and  reported.*5  In  essence,  this  work  Ted  to  the 
selection  of  D-H  Manganin  alloy  and  C-7  epoxy*  as  the  resistive  and  in¬ 
sulator  elements,  respectively,  of  a  high-stress  transducer.  The  response 
characteristics  of  this  transducer  formed  a  major  portion  of  the  ensuing 
investigation.  In  particular,  investigations  of  transducer  pressure- 
resistance  hysteresis,  peak  pressure  calibration,  and  effects  of  insulator 
electrical  and  shock  properties  were  examined.  It  was  found  that  all  these 
could  be  measured  experimentally  with  the  exception  of  insulator  shock 
heating.  Preliminary  calculations  showed  that  C-7  temperatures  may  be 
as  high  as  500°C  at  150  kbar  peak  pressure. 

The  work  reported  in  the  following  sections  is  an  extension  and 
refinement  of  the  work  areas  outlined  above.  One  additional  area  of 
effort  is  described:  matching  the  shock  characteristics  of  the  gage 
to  that  of  proposed  shot  media  by  (1)  decreasing  the  compressibility  of 
the  gage  insulator  and  (2)  imbedding  the  Manganin  directly  in  samples  of 
the  proposed  media. 

• 


Araatreeg  Prefect.,  he,  Venae,  ledieea 


The  purpose  of  this  report  is  primarily  to  describe  the  response 
of  the  Manganin  in  epoxy  transducer  with  respect  to  the  tests  performed 
on  it  and  will  therefore  include  portions  of  work  reported  previously. 

In  addition,  the  most  recent  gage  configurations  and  test  results  are 
described.  Section  II  presents  a  summary  of  work  performed  since  the 
report  of  Reference  4;  Section  III,  a  brief  discussion  of  some  aspects 
of  the  results;  Section  IV  A  describes  the  shock  characteristics  of  C>7; 
and  Appendix  A  and  B,  the  C-7  Hugoniot  temperature  calculations;  IV  B, 
the  shock  response  of  Manganin;  IV  C  the  Hugoniot  matching  work;  IV  D, 
recent  transducer  configurations;  IV  E,  recording  circuits;  and  IV  F,  field 
test  work. 
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II  SUMMARY 


Response  of  the  Manganin  in  epoxy  transducer  has  been  extended  to 
~  290  kbar  peak  pressure  calibration  and  ~  150  kbar  hysteresis  calibra¬ 
tion.  Constant  pressure  (ambient)  measurements  have  been  made  of  the 
temperature  coefficient  up  through  the  melting  temperature,  ~  1020°C. 

The  transformation  from  transducer  profile  recordings  to  profiles  in 
adjacent  test  materials  has  been  examined  by  the  method  of  characteris¬ 
tics.  For  the  numerical  example  chosen,  C-7  epoxy  on  aluminum,  the 
distortion  caused  by  the  presence  of  gage  is  small..  Hugoniot  temperature 
calculations,  completed  for  C-7  epoxy,  shov;  epoxy  temperatures  reaching 
~  2500°K  at  500  kbar.  Wire  temperature  calculations  as  a  function  of 
time  have  not  been  completed. 

The  feasibility  of  making  profile  measurements  directly  in  limestone, 
granite,  and  tuff  with  a  Manganin  wire  has  been  demonstrated.  Complete 
Manganin  response  characteristics  in  these  media  have  not  been  deter¬ 
mined;  however,  measurements  of  peak  pressure  were  about  as  expected. 

An  insulator-Manganin  foil  assembly  has  been  used  successfully  in  an 
aluminum  block  configuration  to  record  peak  pressure  and  profile  for 
~  1  /xsec.  The  mechanism  of  failure  after  this  time  has  not  been 
determined. 


Ill  DISCUSSION 


A.  GENERAL 

A  transducer  capable  of  recording  shock  wave  stress-time  profiles 
must  necessarily  respond  to  very  high  frequency  changes  in  mechanical 
stress.  Consequently,  stress  wavelengths  are  generally -much  less  than 
transducer  dimensions,  and  the  shock  wave  response  of  the  individual  trans¬ 
ducer  components  determines  not  only  the  overall  transducer  response  but 
also  the  actual  profile  being  measured.  Such  is  the  case  with  Manganin 
wire  in  C-7  epoxy,  shown  schematically  in  Fig.  1(a)  and  pictorially  in 
1(b).  The  wire  is  the  transducing  element,  and  the  surrounding  C-7  in¬ 
sulator  acts  as  a  pressure  transformer  between  a  test  specimen  and  the 
wire.  The  stress-time  profile  recorded  depends  on  the  shock  properties 
of  the  insulator  as  well  as  the  piezoresistive  response  of  the  wire. 

Ideally,  the  insulator  should  behave  as  a  fluid  of  low  electrical 
conductivity.  Measurements  of  transducer  response  and  C-7  equation  of 
state  reported  in  the  following  sections  indicate  that,  within  experi¬ 
mental  error,  C-7  behaves  ideally  up  to  at  least  150  kbar.  Measurements 
above  this  pressure  indicate  a  possible  change  in  slope  in  the  shock 
velocity-pressure  velocity  curve.  Frequently  this  is  indicative  of  a 
phase  transformation  of  a  material  and  is  characterized  in  the  pressure¬ 
time  profile  as  a  double  wave.6  However,  a  double-wave  structure  has  not 
been  observed  in  C-7,  although  both  free-surface  velocity7  and  transducer 
records  have  been  recorded  up  to  ~290  kbar. 

Assuming  the  epoxy  does  behave  as  a  fluid,  the  Manganin  wire  peak 
pressure  calibration,  which  has  been  obtained  principally  in  C-7  epoxy, 
must  be  considered  as  applicable  only  when  the  wire  is  under  dynamic 
hydrostatic  pressure.  Should  the  wire  be  used  in  a  material  in  which 
the  stress  is  anisotropic,  the  response  may  be  altered;  static  pressure 
measurements  on  Manganin  have  indicated  this  to  be  true.8  Dynamic  mea¬ 
surements  in  which  the  anisotropy  is  known  have  not  been  performed;  how¬ 
ever,  these  may  now  be  possible  with  the  transducer  configurations  de¬ 
scribed  in  Section  IV  D.  In  fact,  the  universality  of  these  transducers 
may  depend  on  knowing  the  response  of  Manganin  to  anisotropic  stress. 
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FIG.  1  MANGANIN  IN  C-7  EPOXY  TRANSDUCER 


By  referring  the  Manganin  peak  pressure  calibration  to  the  epoxy 
pressure,  there  is  one  additional  restriction:  the  results  may  be  valid 
only  for  the  particular  method  by  which  the  pressure  is  obtained,  i.e., 
a  single  shock  in  C-7.  Experiments  with  Manganin  in  other  insulators, 
notably  doped  epoxy,  have  yielded  calibration  values  within  the  experi¬ 
mental  error  of  those  obtained  with  C-7.  However,  the  compressibilities, 
and  hence  the  paths  by  which  the  final  P-V-E  (pressure- volume-energy ) 
states  are  reached  have  been  very  similar  to  C-7.  Two  possible  excep¬ 
tions  are  the  Hi-D  glass  results,  reported  under  “Peak  Pressure  Calibra¬ 
tion,”  and  Manganin  in  aluminum,  reported  under  "Additional  Gage  Develop¬ 
ment.”  To  date,  sufficient  data  have  not  been  generated  to  distinguish 
between  scatter  and  a  change  in  the  pressure  coefficient  in  either  sys¬ 
tem  when  compared  to  C-7;  however,  this  will  be  studied  further,  both 
theoretically  and  experimentally. 

Values  of  pressure- resistance  hysteresis  have  been  determined  for 
the  combination  of  wire  in  epoxy.  Pressure  cycling  was  performed  by  shock 
compression  and  subsequent  unloading  by  relief  wave  reflection  from  a  free 
surface.  Since  the  expansion  under  these  conditions  is  adiabatic,6  the 
epoxy  and  wire  do  not  return  to  prestressed  temperature  and  volume.  Pos¬ 
sible  contributions  of  the  residual  temperature  to  measured  hysteresis 
values  have  been  assumed  negligible  since  zero-pressure  temperature  co¬ 
efficient  measurements  indicate  insignificant  resistance  change,  <~2%, 
up  to  approximately  the  Manganin  melting  temperature.  To  the  extent  that 
measurements  apply  after  pressure  cycling,  the  quoted  hysteresis  values 
can  be  considered  applicable  to  Manganin. 

B.  TRANSDUCER  APPLICATION 

It  was  mentioned  above  that  the  Manganin  in  C-7  transducer  measures 
a  profile  existing  in  the  C-7  insulator.  Hie  corresponding  peak  pres¬ 
sure  and  pressure  profile  in  an  adjacent  test  material  depend  on  the  shock 
impedances  and  equations  of  state,  respectively,  of  the  two  materials. 

A  transformation  in  peak  pressure  can  be  readily  accomplished  by  the  well 
known  impedance  match  procedure  if  the  appropriate  cross  characteristics 
of  the  materials  are  known.9  A  pressure- time  profile  transformation  is 
somewhat  more  complex,  depending  on  the  quantity  desired;  that  is,  the 
profile  measured  by  the  transducer  is  that  existing  at  the  test  material- 
transducer  interface.  The  profile  in  the  test  material  is  altered  by  the 
waves  reflected  from  the  interface.  In  field  test  work  this  effect  is 
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not  included  in  predicted  profiles.  Hence  an  estimate  of  the  distortion 
caused  by  the  gage  is  required.  This  can  be  accomplished  for  materials 
which  behave  as  fluids  by  a  calculation  using  the  method  of  character¬ 
istics,  10  as  outlined  in  Appendix  A.  A  specific  example  (suitable  for 
comparison  with  future  measured  profiles)  of  a  given  driver  system  was 
chosen  for  analysis.  The  general  procedure  is  discussed  first  step  by 
step  to  demonstrate  the  construction  of  a  characteristic  “net”  from  which 
appropriate  pressures  can  be  assigned  to"  regions  of  an  adj acent- test 
material  as  represented  in  distance-time  (x,t)  space. 

For  the  materials  used,  the  net  effect  of  interaction  of  relief  waves 
is  a  decrease  in  the  rate  of  decay  of  pressure  from  that  which  would  ex¬ 
ist  if  the  gage  material  were  identical  to  the  test  material.  In  the  ex¬ 
ample  chosen,  a  true  rate  of  decay  (defined  as  that  existing  in  an  all 
aluminum  system  or  if  the  gage  insulator  were  matched  to  the  aluminum) 
can  be  estimated  by  extending  the  unrefracted  relief  waves  to  the  right 
in  Fig.  A-l  and  through  the  plane  of  the  (virtual)  interface.  Pressure 
vs.  time  at  the  epoxy-aluminum  interface  and  at  the  virtual  interface, 
as  determined  from  gage  data,  are  given  in  Table  I.  Here  the  time  incre¬ 
ments  in  column  2  result  from  the  arbitrarily  chosen  time  increments, 

Table  I 

PRESSURE-TIME  PROFILE  DISTORTION  RESULTS 


PRESSURE 

(kb«r) 

TIRE  OF  RELIEF  WAVE  ARRIVAL  <m*«c)  | 

Aliaiaia, 

C-7  Iittrfica 

148 

0 

0 

140 

0.44 

0.42 

129 

0.84 

117 

1.30 

1.23 

106 

1.61 

1.53 

% 

2.00 

1.85 

90 

2.39 

2.30 

79 

2.77 

2.69 

69 

3.16 

3.01 

0.4  ^sec,  used  in  the  analysis  of  Appendix  A.  Zero  time  refers  to  shock 
arrival  at  the  aluminum  C-7  interface.  Hie  time  intervals  shown  in  col¬ 
umn  3  are  those  obtained  graphically  without  the  gage  present,  i.e.. 
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in  an  all  aluminum  medium.  It  can  be  seen  that  the  difference  is  less 
than  the  desired  iO%  accuracy  set  as  a  gage  goal.  For  this  particular 
example,  distortion  in  the  profile  caused  by  the  presence  of  the  gage 
would  be  small.  For  estimates  of  other  materials,  the  Hugoniots  and  re¬ 
lief  adiabats,  from  which  the  sound  velocity  can  be  computed,  must  be 
known . 
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IV  EXPERIMENTAL  PROCEDURE  AND  RESULTS 

A.  SHOCK  CHARACTERISTICS  OF  C-7  EPOXY 

Since  the  integrity  of  the  Manganin  in  C-7  transducer  depends  greatly 
on  the  characteristics  of  the  C-7»  it  is  necessary  that  these  be  examined 
in  considerable  detail.  The  procedure  followed  in  the  present  work  has 
been  first  to  measure  the  C-7  Hugoniot  equation  of  state;  then,  by  using 
these  data,  to  calibrate  the  Manganin  wire  for  peak  pressure  response, 
hysteresis,  etc.  In  addition,  since  shock  heating  accompanies  shock  com¬ 
pression,  a  program  was  initiated  to  calculate  Hugoniot  temperatures  and 
possible  effects  on  wire  resistance.  Temperature  calculations  at  pressures 
up  to  500  kbar  have  been  completed;  calculations  of  wire  temperature  as  a 
function  of  time  have  not  been  started. 

1.  Hugoniot  Measurements 

The  Hugoniot  of  C-7  has  been  reported  previously4  for  pressures  up 
to  150  kbar.  This  has  been  extended  to  220  kbar  by  conventional  “in 
contact"  explosives  and  to  290  kbar  by  flying  plate  techniques.  11  A  com¬ 
plete  summary  of  Hugoniot  data  is  presented  in  Table  II.  The  data  up  to 
~  170  kbar  can  be  represented  by  an  equation  of  the  form: 

P  *  0.07816 m  +  0.1956m2  +  0. 2214m3  (mbar) 

where 

M  =  —  "  1,  Mn  ■  1.19  g/cc  .  (1) 

A  fit  in  the  range  from  170  to  the  highest  pressure  obtained,  370  kbar, 
has  not  been  made  because  of  insufficient  data. 

2.  Hugoniot  Temperature  Calculations 

Shock  loading  of  the  Manganin  wire  imbedded  in  expoxy  causes  a  change 
in  wire  temperature  due  to  direct  shock  heating  and  conduction  heating  from 
the  surrounding  epoxy.  Since  a  pressure  profile  and  not  a  pressure- 
temperature  profile  is  desired,  an  estimate  of  epoxy  temperature  and  its 
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Table  II 

SUMMARY  OF  C-7  HUGONIOT  DATA 


PRESSURE 

V, 

V . 

(kbit) 

(■■/fisec ) 

P 

(a«//xaec) 

METHOD 

37.7 

3.86 

0.83 

F. S.  *  wedge 

40.9 

3.97 

45.9 

-  4.06 

51.0 

4.15 

55.5 

4.24 

1.10 

60.3 

4.34 

1.17 

65.4 

4.42 

1.24 

70.4 

4.51 

1.31 

75.4 

4.60 

1.38 

81.0 

4.60 

1.45 

57 

4.38 

1.12 

Imp.  match  (Al) 

4.33 

145 

5.76 

2.J4 

Imp.  match, (Al) 

F.S. ,  F.P.* 

274.5 

7.19 

3.24 

74 

4.70 

1.33 

F.  S. 

74 

4.70 

1.33 

Imp.  match  (Mg) 

170 

6.2 

2.34 

Imp.  match  (Mg) 

170 

6.2 

2.34 

F.S. 

175 

190 

6.34 

6.34 

2.39 

2.57 

Img.  match  (Mg) 

237 

6.79 

2.99 

F.S. 

225 

6.79 

2.88 

Imp.  match  (Mg) 

F7b.,  F.S. 

370 

7.88 

3.98 

368 

7.88 

3.98 

F. P. ,  Imp.  match  (Al) 

293 

8.06 

3.09 

F.P.,  F.S. 

259 

6.46 

3.39 

F.P.,  F-S. 

20.7 

0.49 

Gage,  ■•aiming 

k  «  0. 29%/kbar 

Fra*  aurfaca  velocity. 
*  Flyiaf  plate. 


effect  on  the  wire  is  necessary.  Calculations  of  the  temperature- time 
profile  in  the  wire  are  being  made  in  two  steps;  calculation  of  epoxy 
Hugoniot  temperatures  and  conduction  heating  of  the  wire  by  the  epoxy. 

The  first  has  been  completed  and  is  outlined  in  Appendix  B.  The  calcu¬ 
lation  is  made  by  assuming  that  the  epoxy  equation  of  state  can  be  repre¬ 
sented  by  the  Mie-Grflnei sen  equation  with  the  zero  pressure  Grfineisen 
constant,  T,  calculated  from  elastic  and  thermal  properties.  At  elevated 
pressures,  values  of  T  are  used  which  are  obtained  from  one  measurement 
and  its  estimated  limits  of  error.  (See  Appendix  B. )  An  isentrope  centered 
at  the  foot  of  the  Hugoniot  is  calculated  by  an  iterative  process  from  the 
measured  Hugoniot  and  several  values  of  F\  Hugoniot  temperatures  are  then 
generated  by  using  values  of  specific  heat  at  constant  volume  which  were 
calculated  from  specific  heat  measured  isobarically  as  a  function  of 
temperature.  The  applicability  of  the  Mie-Gr&neisen  equation  of  state  to 
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C-7  plastic  and  the  dependence  of  the  specific  heat  on  volume  are  not 
known.  However,  the  calculations  were  made  to  obtain  only  an  estimate  of 
plastic  shock  temperatures  from  which  wire  temperature  could  be  calculated 
and  are  felt  to  be  sufficient  for  this  purpose. 

Results  of  the  Hugoniot  temperature  calculations  are  shown  in  Fig.  B-5. 
The  significance  of  the  various  values  shown  for  T  and  the  regions  labeled 
I  and  II  are  explained  in  Appendix  B.  It  can  be  seen  that  epoxy  tempera¬ 
tures  are  apparently  high  even  at  modest  pressures,  being 1000 °K  at 
200  kbar. 

3.  Electrical  Properties 

In  general,  the  conductivity  and  charge  generation  measurement  tech¬ 
niques  outlined  in  Reference  4  have  been  applied  to  the  extended  pressure 
range  up  to  29p  kbar.  An  exception  is  a  parallel  plate  method  in  which  a 
more  quantitative  value  of  conductivity  was  obtained.  The  configuration 
is  shown  in  Fig.  2(a). 

Any  polarization  disturbances  induced  on  one  plate  must  likewise  be 
induced  on  the  other  because  the  plates  are  identically  positioned  with 
respect  to  the  shock  front.  The  external  circuit  provides  nearly  identical 
paths  to  ground  from  each  plate,  and  consequently  polarisation-induced 
signals  in  the  direction  of  shock  propagation  will  not  be  present  in  Vl2> 
the  voltage  displayed  on  the  scope.  F12  is  proportional  to  the  current 
which  flows  through  the  sample  because  of  the  applied  voltage.  This  yields 
the  resistance  of  the  sample  from  which  the  resistivity  may  be  easily 
computed. 

Another  pair  of  plates  is  placed  near  the  first  pair  and  is  connected 
to  a  similar  external  circuit;  however,  no  voltage  is  applied  to  the  sample. 
This  circuit  is  monitored  as  a  check  on  the  polarization  signals  produced 
normally  to  the  direction  of  shock  propagation. 

An  experiment  with  this  arrangement  has  been  performed  using  C-7  epoxy 
plastic  as  the  insulator  and  aluminum  foil  plates  5  mm  by  15  mm  with  a 
separation  of  1  mm.  This  gave  a  resistance  of  1.9  *  10s  Q  at  140  kbar 
immediately  after  shock  front  passage  and  corresponds  to  a  resistivity  of 
1.6  x  10“  /iflcm,  neglecting  fringing  effects  at  the  edges  of  the  plates. 
(Inclusion  of  a  correction  for  fringe  effects  would  raise  this  figure 
slightly.)  Thus  the  quoted  value  is  a  minimum  resistivity  which,  when  com¬ 
pared  to  the  resistivity  of  Manganin,  ~  48  /xfi  cm,  is  essentially  infinite. 
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(b)  tA-im-ut 

FIG.  2  INSULATOR  RESISTIVITY  MEASUREMENT 

(a)  Elactrods  Configuration 

(b)  Circuit  Sc  Somatic 

Induced  currents  in  the  Manganin  wire,  attributed  to  epoxy  dipole  re- 
oreintation  under  shock  loading,  have  been  observed  in  transducer  records 
during  shock  transit  through  the  epoxy  prior  to  reaching  the  wire.  Hauver2 
has  investigated  this  phenomenon  extensively  in  polar  plastics  and  has 
generated  quantitative  data  relating  induced  currents  to  peak  pressure. 

In  the  present  work  induced  currents  are  <  4  percent  of  the  measuring 
current  at  the  highest  calibration  pressure  obtained,  290  ±  30  kbar. 

Induced  currents,  while  small,  do  create  some  difficulty  in  the  mea¬ 
surement  of  Manganin  resistance  change.  The  measurement  is  made  with  a 
constant  current  power  supply  system  (see  Section  IV  E);  however,  induced 
currents  cause  a  change  in  current  just  prior  to  shock  arrival.  This 
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appears  on  transducer  records  as  a  small  step  preceding  the  pressure  wave. 

In  the  present  work  it  has  been  assumed  that  the  current  at  the  instant  of 
shock  arrival  is  the  sum  of  power  supply  and  induced  currents.  Any  error 
introduced  by  this  assumption  is  relatively  small,  as  can  be  seen  from  the 
order  of  magnitude  of  the  current  as  quoted  above. 

t 

B.  SHOCK  RESPONSE  OF  MANGANIN 

1.  MANGANIN-HUGONIOT  EQUATION  OF  Sf ATE 

The  Hugoniot  equation  of  state  of  Manganin  has  been  measured  pre¬ 
viously4  up  to  359  kbar  peak  pressure.  The  results  are  included  here 
(Table  III)  for  comparison  with  the  piezoresistance  data  presented  in 
Table  IV  and  curve  of  Fig.  3.  It  can  be  seen  that  there  are  no  indica¬ 
tions  of  a  phase  change  (characterized  by  an  inflection  point  in  the 
Hugoniot6).  A  phase  change  would  raise  the  possibility  of  a  discontinuity 
in  the  pressure- resistance  relationship.  Peak  pressure  calibration  reported 
below  has  shown  no  discontinuity  over  essentially  the  same  pressure  range. 

2.  Peak  Pressure  Calibration 

The  resistance  of  Manganin  as  a  function  of  peak  pressure  has  been 
measured  principally  for  one  alloy,  Driver-Harris  (nominal  composition, 

84%  Cu,  12%  Mn,  4%  Ni ) .  A  summary  of  similar  work  on  other  Manganin  alloys 
has  been  presented  elsewhere. 12  The  results  discussed  below  pertain  only 
to  the  Driver-Harris  alloy. 

In  a  pressure  range  from  20  to  200  kbar,  "in  contact"  explosive 
driver  systems  have  been  used  to  generate  dynamic  pressures.  A  typical 
assembly  of  this  type  is  shown  in  Fig.  4.  A  range  of  peak  wire  pressures 
was  obtained  by  using  insulator  and  driver  plate  materials  of  different 
shock  impedances,  by  using  high  explosives  of  different  Chapman- Jouget 
pressures,  and  by  varying  the  explosive  pad  thickness  and  hence  the  peak 
pressure  attenuation  due  to  the  Taylor  relief  wave.  The  majority  of  work 
has  been  done  with  one  insulator,  C-7  epoxy,  since  the  change  in  resistance 
of  Manganin  with  pressure  has  not  yet  been  shown  to  be  independent  of  the 
internal  energy  state  of  the  Manganin  at  a  given  pressure.  For  this  reason, 
peak  pressure  vs.  resistance  values  are  considered  uniquely  determined  by 
the  method  by  which  the  peak  pressure  was  obtained;  for  example,  direct 
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Table  III 

HUGONIOT  DATA,  MANGANIN  ALLOY, 
P0  =  8.46  g/cc 


PRESSURE 

(kbar) 

SHOCK 

VELOCITY 

(aa/paee ) 

PARTICLE 

VELOCITY 

(an/psec) 

RELATIVE 

VOLUME 

1  -  (v'v0) 

METHOD  OF 
OBSERVATION 

56 

4.02 

0.165 

0.041 

Free  surface 
(inclined  mirror) 

107 

4.36 

0.294 

0.067 

Impedance  match 

157 

4.52 

0.413 

0.091 

Free  surface 
(inclined  mirror) 

254 

4.88 

0.616 

0.126 

Free  surface 
(inclined  mirror) 

346 

5.15 

0.780 

0.150 

Free  surface 
(inclined  mirror) 

360 

5.18 

0.820 

0.159 

Free  surface 
(inclined  mirror) 

Table  IV 

PRESSURE-RESISTANCE  DATA 
MANGANIN  WIRE 


PRESSURE 

(kbar) 

INSULATOR 

20  *  1 

24  4  2 

07 

59  t  2 

74  i  2 

134  t  4 

166  1  $ 

169  t  5 

17$  t  $ 

212  t  S 

220  t  6 

120  1  $ 

281  i  10 

360  t  30 

135  1  S 

Hi-D 

glass 

279  t  10 

129  i  8 

Alaaiaam* 

290  ±  30 

07 

DIRECT  (A) 

OR  REFLECTED  SHOT 
(A)  SHOCK  NO. 


£JVI>0  (*) 


Oacilloicope  No. 


iJVH  q  1  or 

X 


PRESSURE  SOURCE 


Light  gaa  gas 

“la  coatact*  txploaivt 


16.9 
20.6 
3B.6 
47.0 
47.0 
48.1 

52.9  ±  1.7 
$5.3  i  1.9 


61.  S 

33.2  1  1.4 


69.0  I  Flying  plate 


H-im-MC 


FIG.  4  SCHEMATIC  CONFIGURATION  OF  DRIVER  SYSTEM, 
TRANSDUCER,  AND  RECORDING  CIRCUITS 


shock  vs.  reflected  shock,  and  also  according  to  the  medium  in  which  the 
wire  was  imbedded.  The  results  of  the  most  recent  calibration  shots  are 
presented  in  Table  IV  and  the  curve  of  Fig.  5. 

Calibration  pressures  have  been  determined  by  measuring  insulator 
pressures  by  either  the  impedance  match  procedure  or  a  free  surf ace- shock 
velocity  measurement.^  The  Manganin  wire  is  placed  as  close  to  the  plane 
of  free-surface  velocity  measurement  or  to  the  driver  plate-insulator  in¬ 
terface  as  possible,  usually  within  1  mm.  In  this  way,  errors  introduced 
by  peak  pressure  attenuation  are  minimized. 

The  calibration  of  Manganin  in  C-7  as  a  transducer  is  statistical  in 
that  the  transducer  is,  of  course,  destroyed  in  use.  A  brief  discussion 
of  the  reliability  of  the  calibration  as  well  as  accuracy  is  therefore 
appropriate. 

In  general,  it  was  found  that  if  several  wires  were  imbedded  at  the 
same  plane  in  C-7  epoxy,  the  change  in  resistance  agreed  to  within  the 
reading  accuracy  of  the  recording  oscilloscopes.  Under  the  best  condi¬ 
tions;  that  is,  with  large  percentage  changes  in  resistance,  this  accuracy 
was  ±1.5  percent,  deteriorating  to  ±5  percent  at  the  lowest  resistance 
change.  The  accuracy  of  the  corresponding  pressure  measurement  varied 
inversely  with  pressure,  being  typically  ±5  percent  at  the  low  pressures 
and  as  high  as  ±10  percent  at  the  highest  direct  shock  pressure  reached. 

The  pressure  points  above  220  kbar  listed  in  Table  IV  were  obtained, 
as  indicated,  by  one  or  a  combination  of  several  techniques;  flying  plates 
reflected  shock  wave,  and  Hi -D  glass,  or  aluminum  media.  The  technique 
used  in  the  aluminum  shots  will  be  discussed  in  Section  IV  D.  In  the  case 
of  the  flying  plate  shots,  the  tabulated  result  represents  one  successful 
shot  out  of  three  attempts.  In  this  one  shot,  although  two  separate  trans 
ducers  gave  resistance  changes  within  3  percent,  optical  measurements  of 
free-surface  and  shock  velocity  yielded  the  large  uncertainty  shown.  This 
ambiguity  in  optical  measurements  resulted  from  the  non-planar  shock 
arrival.  The  pressure  error  limits  shown  are  based  on  errors  incurred  in 
data  reduction. 

The  reflected  wave  shot  yielded  more  consistent  results,  although  a 
comparison  with  direct  shock  values  can  not  be  made  at  this  time,  because 
of  the  uncertainty  in  the  PVE  states,  as  described  above.  The  configura¬ 
tion  of  these  shots  is  shown  in  Fig.  6.  The  peak  pressure  of  the  brass 
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FIG.  5  CHANGE  IN  RESISTANCE  AS  A  FUNCTION  OF  PRESSURE,  D-H  MANGANIN  WIRE 


backing-plate  was  obtained  by  a  free-surface  velocity  measurement  from 
which  the  C-7  pressure  was  obtained  by  an  impedance  match  procedure.  At¬ 
tenuation  in  the  reflected  and  incident  wave,  both  in  the  C-7  and  in  the 
thin  brass  plate,  has  been  neglected  for  two  reasons:  the  close  proximity 
of  the  wire  to  the  epoxy-brass  interface  and  the  slow  rate  of  change  -of 
pressure  vs.  time  as  indicated  by  the  wire  record. 

It  can  be  seen  from  the  curve  of  Fig.  5  that,  for  the  direct  shock 
epoxy  points,  there  appears  to  be  a  change  in  slope  occurring  in  the  region 
of  150  kbar.  the  reflected  wave  and  Hi-D  data,  although  limited,  do  not 
show  this.  It  is  suspected  that  the  change  in  slope  may  be  due  to  the 
epoxy  and  not  the  Manganin;  however,  C-7  Hugoniot  results  do  not  reflect 
this.  Considerably  more  data  and  more  accurate  pressure  measurements  are 
required  to  resolve  the  pressure- resistance  dependence  of  Manganin  further. 
Within  the  accuracy  desired  for  field  gages,  the  change  in  resistance  of 
Manganin  with  peak  pressure  in  C-7  epoxy  and  under  direct  shock  loading 
can  be.  represented  by: 

=  iltj  •  P(kbar)  for  P  <  ~  150  kbar  (2) 

-  0.16  +  k2  •  P(kbar)  for  150  <  P  <  290  kbar  (3) 

■  0.29  ±  0.02  (%/kbar),  Jk,  -  0.18  ±  0.01  (%/kbar)  . 

The  nominal  value  of  fcj  is  about  4  percent  less  than  that  previously 
quoted.  4,13  The  principal  reason  for  this  is  a  change  in  the  recording 
system  (explained  in  Section  IV  E)  to  permit  more  accurate  results  and 
higher  frequency  response. 

3.  Temperature  Dependence  of  the  Pressure  Coefficient 

The  entent  to  which  wire  heating,  either  directly  by  shock  or  through 
conduction  from  the  insulator,  affects  the  measured  resistance  change  is 
difficult  to  determine.  Static  pressure  measurements13  have  indicated  that 
the  effect  is  negligible.  Dynamic  measurements,  reported  previously,4  have 
been  indirect  in  that  wire  heating  was  accomplished  before  shock  compression, 
whereas  to  simulate  actual  conditions  in  an  insulator,  heating  should  be 
performed  after  shock  compression.  The  results  reported  in  Reference  4 
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show  that  up  to  500°C  and  145  kbar  the  pressure  coefficient  is  indepen* 
dent  of  temperature  within  experimental  error.  Recent  work  has  not 
extended  this  measurement  but  has  been  confined  to  measuring  the  ambient 
pressure  temperature  coefficient  up  through  the  melting  point  of  Manganin 
(~1000°C),  thereby  obtaining  at  least  an  indication  of  its  behavior. 

The  results  are  shown  in  the  resistance  vs.  temperature  curve  of  Fig.  7. 

The  largest  change  in  resistance  observed  up  to  the  melting  point  was 
~~2  percent.  The  change  in  resistance  between  500°C,  the  temperature 
reached  in  the  above  referenced  work,  and  melting  temperature  is  negligible 
Consequently,  on  the  basis  of  the  ambient  pressure  temperature  coefficient, 
there  is  no  indication  of  a  change  in  the  pressure  coefficient  with  tem¬ 
perature,  i.e.,  in  shock- loading  Manganin,  no  change  in  the  coefficient 
due  to  an  increase  in  temperature  by  shock  heating  is  expected.  Similarily 
if  the  ambient  pressure  temperature  coefficient  can  be  applied  to  elevated 
pressure  results,  it  would  be  expected  that  conduction  heating  of  the  wire 
would  not  alter  the  pressure-induced  resistance-time  profile  until  perhaps 


FIG.  7  RESISTANCE  AS  A  FUNCTION  OF  TEMPERATURE, 
D-H  MANGANIN  ALLOY 
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the  melting  temperature  of  Manganin  were  reached.  Since  this  is  likely 
to  be  well  above  the  zero-pressure  value14  and  is  difficult  to  predict, 
and  since  conduction  heating  has  not  yet  been  calculated,  the  actual 
effect  of  temperature  can  only  be  inferred  from  the  above  mentioned  results 
and  measured  profiles  reported  below. 

4.  Pressure-Resistance  Hystersis 

Static  pressure  work  with  Manganin  has  shown  the  pressure  coefficient 
to  decrease  slightly  with  repeated  pressure  cycling  until  a  stable  value 
is  reached.  15  In  addition,  the  static  pressure  coefficient  measured  for 
increasing  pressure  is  within  a  few  percent  of  that  measured  for  decreasing 
pressure. 14  These  static  measurements  were  made  isothermally  although  sepa¬ 
rate  temperature  cycling  has  shown  a  stablizing  effect  on  the  pressure 
coefficient.  **  In  shock  wave  studies  it  is  more  difficult  to  isolate  and  to 
measure  a  purely  pressure-resistance  hysteresis  effect  since  wire  compres¬ 
sion  is  accomplished  by  shock  heating  and  expansion  occurs  adiabatically. 
Moreover,  measurements  performed  in  an  insulating  medium  are  influenced 
by  the  thermal  properties  of  the  insulator. 

The  hysteresis  results  reported  here  have  been  made  with  respect  to 
the  Manganin  and  C-7  epoxy  insulator  as  a  combination.  Two  principal 
methods  have  been  used:  1)  subjecting  wire  and  epoxy  to  a  sudden  known 
decrease  in  pressure  shortly  after  peak  pressure,  and  2)  comparing  the 
pressure- time  curve  recorded  in  test  shots  with  one  calculated  from 
pressure-distance  data  measured  on  the  same  shot.  It  is  thought  that  in 
the  case  of  a  step  decrease,  the  wire  expansion  is  adiabatic,  i.e.,  the 
wire  is  insulated  from  the  epoxy,  whereas  in  the  slowly  decaying  wave  heat 
is  conducted  to  the  wire  from  the  epoxy.  However,  conduction  calculations 
have  not  yet  been  made.  In  both  cases,  described  below,  measurements, 
although  few,  have  been  sufficient  to  establish  the  magnitude  of  the 
hysteresis  effect. 

The  experimental  system  for  generating  a  step  decrease  in  pressure 
is  shown  in  Fig.  8.  It  is  similar  to  the  reflected  wire  shots  except  that 
materials  of  lower  shock  impedance  were  placed  immediately  behind  the  sen¬ 
sitive  length  of  wire  (between  voltage  taps,  see  Section  IV  E). 

Two  series  of  shots  have  been  conducted:  pressure  cycling  Manganin 
in  C-7  to  various  peak  pressures  and  subsequent  relief  to  atmospheric,  and 
pressure-cycling  Manganin  in  doped  epoxy  insulators  (see  Section  IV  C)  to 
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FIG.  8  TRANSDUCER  CONFIGURATION, 
RELIEF  WAVE  SHOTS 


various  peak  pressures  and  relieving  to  intermediate  pressures  by  means 
of  C-7  epoxy  backing.  The  results  are  shown  in  Tables  V  and  VI.  The  in¬ 
termediate  pressure  results  are  less  accurate  than  those  at  zero  pressure 
for  two  reasons:  peak  pressure 
and  relief  pressure  values 
are  subject  to  the  rather 
large  uncertainties  in  the 
Hugoniots  and  relief  adiabats 
of  the  doped  epoxies,  and 
quoted  values  of  measured 
decreasing  pressure  change 
were  obtained  by  using  a  pres¬ 
sure  coefficient  obtained 


Table  V 

RELIEF  WAVE  CALIBRATION  DATA, 
INTEMEDIATE  PRESSURES 


INSULATOR 

MEASURED 
Oft/R  (*) 

MEASURED 

-HP 

(kkir ) 

EXPECTED 

-OP 

(kbar) 

1 

2 

67*  Si,  33*  C-T 

C-7 

7.3 

23 

26 

67*  Hi-D,  33*  07 

C-7 

S.S 

19 

17 

77*  Hi-D,  23*  C-7 

C-7 

13.1 

42 

47 

77*  Hi-D,  23*  C-7 

C-7 

17.4 

S6 

SO 
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Table  VI 

RELIEF  WAVE  CALIBRATION  DATA, 
ZERO  PRESSURE 


for  increasing  pressure  change.  It  has  not 
yet  been  shown  that  the  two  are  equivalent. 
However,  the  intermediate  pressure  results 
do  indicate  that  hysteresis  is  insignificant 
compared  to  a  10  percent  accuracy  goal./ 

The  zero-pressure  relief  wave  Shots  show 
a  definite  change  in  zero-pressure  resistance 
after  pressure  cycling.  This  change  appears 
to  be  a  function  of  peak  pressure  as  shown 
in  Fig.  9.  The  changes  in  resistance  have 
not  been  corrected  for  possible  residual 
temperature  effects  since  measurement  of  the  temperature  coefficient  of 
resistance,  outlined  above,  indipates  the  correction  to  be  negligible. 

The  magnitude  of  the  residual  change  after  cycling  to  150  kbar  approaches 
the  desired  overall  accuracy  of  the  transducer  system.  Additional  work 
of  this  type  is  required  to  obtain  more  quantitative  results  and  to  extend 
the  range  of  measurement. 


FIG.  9  PRESSURE-RESISTANCE  HYSTERESIS  OF  MANGANIN 
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CURRENT  AND  VOLTAGE  LEADS 


FIG.  10  CONFIGURATION  OF  P(t),  P(x)  SHOTS 


Response  to  a  gradually  decaying  wave  has  been  examined  by  the  system 
shown  in  Fig.  10.  By  using  several  wires  spaced  a  known  distance  apart  in 
the  direction  of  shock  propagation,  the  measured  pressure-time  profile  at 
any  wire  can  be  compared  with  the  profile  calculated  from  peak  pressure 
decay  with  distance  as  measured  by  subsequent  wires. 


This  can  be  clarified  by  reference  to  Fig.  11,  distance  vs.  time, 
where  at  time  t0,  the  shock  front  is  at  plane  *0  and  peak  pressure  is  Pfl. 
At  time  1 1 ,  the  front  has  moved  to  »1P  and  the  peak  pressure  is  now  P. 
There  is  a  distance  behind  the  shock  front  where  the  relief  wave  pressure 
is  P.  If  we  call  this  distance  coordinate  and  corresponding  time  co¬ 
ordinate  t  we  can  solve  for  t  as  follows: 
p  p 


t 


l 


u  +  a 

p 


(4) 
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and 


UP  o 

where  up  is  the  mass  velocity,  a  is  the  local  sound  velocity  at  pressure 
P,  and  up0  is  the  average  mass  velocity  between  x0  and  xp . 

Since 

*1  ~  *0  =  "  ‘o*  (6) 

where  U  is  the  average  shock  velocity  between  *j  and  *0,  adding  (4)  and 
(5)  and  equating  to  the  right  hand  side  of  (6)  gives 


(Ut  ~  up  ~  a) 
l  ~ 

(uPo~up  "  •> 


(7) 


Subsequent  wires  measure  a  peak  pressure  and  a  corresponding  average 
shock  velocity,  Ut .  Values  of  up0  are  known  from  Hugoniot  data.  The 


FIG.  11  DISTANCE-TIME  DIAGRAM, 
RELIEF  WAVE  CALIBRATION 
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sound  speed,  a,  -can  be  calculated  from  the  difference  in  slope  between  the 
Hugoniot  and  adiabat  at  a  given  pressure,  P,  and  volume,  V,  as  given  by 6 

a2  =  c2(l  -  A)  ,  (8) 

where  a  is  the  adiabatic  sound  speed,  c2  =  dP/dp  along  the  Hugoniot  and  A 
i s  defined  by : 6 

A  =  7  P<-V o  ~  V)[l  ‘  <0  ‘  «,)2/c2 1  (9) 

r  =  Griineisen  ratio  =  V/cu(3p/3T)  B,  cv  is  the. specific 
heat  at  constant  volume. 

V  =  specific  volume  55  1/p  cc/g 

U  *  shock  velocity  at  P,  V 
Up  ■  mass  velocity  at  P,  V. 

The  ratio  oT_c/d  as  a  function  of  pressure  for  C-7  epoxy  is  shown 
in  Fig.  12. 

Three  shots  were  made  with  a  gage  configuration  as  shown  in  Fig.  10, 
two  at  '''ISO  kbar  peak  pressure  in  the  first  wire  and  one  at  ''‘63  kbar  peak 
pressure.  In  all  shots  the  difference  between  measured  and  calculated 
profiles  for  each  wire  was  not  greater  than  10  percent.  A  tabulation  of 
predicted  and  measured  values  for  one  of  the  high  pressure  shots  is  shown 
in  Table  VII,  where  Pn  f  is  the  predicted  pressure  at  wire,  n,  from  peak 
pressure  measurements  at  wire,  r,  etc.,  and  Pm  is  the  meecured  pressure 
at  time,  t  +  *,  determined  by  the  procedure  outlined  above.  The  longest 
interval  available  for  comparison  is  under  2  psec.  Longer  times  require 

Table  Vll 

PC t )  DATA  DETERMINED  FROM  P(x)  MEASUREMENTS, 

EIGHT- INCH  PLANE  WAVE  LENS  (P  •  80)  PLUS 
1-INCH  OOMPBPAD 
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a  prohibitively  large  plane-wave  explosive  system  to  prevent  lateral  re¬ 
lief  waves  from  distorting  the  profile.  However,  the  measured  response 
is  within  the  desired  accuracy. 

C.  HUGONIOT  MATCHING 

1.  General 

Although  the  majority  of  equation  of  State  measurements  have  been  on 
C-7  epoxy,  Hugoniot  points  have  been  obtained  on  a  number  of  other  less 
compressible  insulators.  The  goal  of  this  work  has  been  to  provide  a 
spectrum  of  transducer  insulators  of  shock  characteristics  closely  matched 
to  expected  field  test  environments.  The  initial  work4  was  directed  toward 
matching  shock  impedances  to  reduce  peak  pressure  and  pressure  profile  errors 
introduced  by  uncertainties  in  the  shock  characteristics  of  proposed  test 
media.  The  work  reported  below  has  been  directed  toward  matching  transducer 
insulator  compressibilities  to  those  of  its  surroundings.  By  doing  this, 
it  is  expected  that  distortionless  recording  for  long  durations  can  be 
achieved  since  disturbances  originating  at  gage  lateral  boundaries  will 
be  eliminated.  Two  methods  have  been  pursued:  adding  less  compressible 
materials  to  expoxies  ("doped'’  epoxies)  and  constructing  transducers  from 
native  field  shot  media. 

2.  Hugoniot  Matching 

To  accomplish  Hugoniot  matching,  a  method  of  predicting  compressibil¬ 
ities  of  mixtures  as  a  function  of  pressure  is  first  required.  From  this, 
it  is  possible  to  obtain  an  expression  for  the  compressibility-pressure 
relationship  required  of  an  additive  to  a  given  epoxy  in  order  to  obtain  a 
given  Hugoniot.  Since  considerable  Hugoniot  data  had  been  generated  on 
doped  epoxies, 4  the  data  were  fitted  to  existing  synthesis  methods  to  de¬ 
termine  their  applicability  to  mixtures.  This  portion  of  the  work  has 
been  completed  and  is  presented  below.  The  results  have  not  been  applied 
to  predicting  the  compressibility  characteristics  of  materials  required  to 
match  shot  media. 

The  method  of  Lyakhov17  as  described  by  Chabai  18  appears  to  work 
reasonably  well  when  applied  to  Hi-D  in  C-7  and  W  in  C-7.  As  described 
in  Reference  17,  the  mixture  compressibility  at  a  given  pressure  is  as¬ 
sumed  equal  to  the  sum  of  the  constituent  compressibilities.  A  weighting 
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factor  for  each  constituent  is  used  according  to  constituent  volume  in 
the  mixture,  i.e. , 


N 

T)(P)  =  2  K.  •  7J.(P)  (10) 

«  =  1  1  * 

where  P  is  shock  pressure,  tj  =  p/pQ,  the  ratio  of  shocked  to  unshocked 
density,  and  K  is  the  volume  weighting  factor.  For  the  two-component 
epoxy  system,  the  mixture  compressibility  7j^(P)  equals 

VP)  =  «VP>  +  (1  "  *)VP)  (11) 

where  M,  E,  and  A  refer  to  mixture,  epoxy,  and  additive  material, 
respectively. 

The  results  obtained  with  W  in  C-7,  Hi-D  in  C-7,  and  Si02  in  C-7  are 
shown  in  Table  VIII.  The  calculated  shock  velocities  of  the  mixture  are 
obtained  from  the  Bankine-Hugoniot  conservation  equations 19  and  in  this 
case  are  given  by 


V,  =  /t,(P  -  P0)/(p  -  p„)  (12) 

The  summation  appears  to  yield  good  correlation  between  Ut  and  PQV^,  cal¬ 
culated  and  observed,  for  W  in  C-7  and  Hi-D  in  C-7.  In  the  case  of  Si02 
in  C-7,  the  equation  of  state  of  fused  quartz  was  assumed  to  apply  to  the 
Si02  actually  used.  The  correlation  was  fair  but  not  so  good  as  for  the 
other  two  cases. 

3.  Limestone,  Tuff,  and  Granite  Transducers 

Construction  of  transducers  composed  of  native  shot  media  has  been 
accomplished  for  three  materials,  limestone,  tuff,  and  granite.  The 
general  configuration  is  shown  in  fig.  13.  It  consists  of  a  block  as¬ 
sembly  cemented  together  with  C-7  epoxy.  A  conventional  four-lead 
(voltage  and  current),  3-mil-diameter  Manganin  wire  assembly  lies  on  the 
block  marked  No.  1  and  under  block  No.  4.  The  voltage  leads,  Cu  foils 
0.004  inch  thick  *  0.125  inch  wide,  are  brought  to  the  cable  end  of  the 
gage  (top)  between  blocks  Nos.  4  and  2  and  similar  current  leads  between 
blocks  Nos.  3  and  4,  and  blocks  Nos.  5  and  4-  Shown  also  is  a  conductivity 
test  foil  extending  from  the  plane  of  the  wire  to  the  cable  end  of  the  gage 
(between  block  Nos.  5  and  6). 


34 


BLOCK  DIMENSIONS  —  inches 
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FIG.  13  LIMESTONE  BLOCK  GAGE  ASSEMBLY,  SCHEMATIC  VIEW 


A  pressure- time  record  from  one  limestone*  shot  using  this  assembly 
is  shown  in  Fig.  14.  The  upper  trace  is  that  of  the  conductivity  test 
probe  with  a  2-fi  external  shunt  resistor  and  associated  constant  current 
power  supply.  It  can  be  seen  that  the  limestone  remains  an  insulator  with 
respect  to  the  2-fl  shunt.  Hie  signal  on  the  bottom  trace  is,  therefore-, 
that  of  the  Manganin  wire  which  is  a  direct  measure  of  the  pressure-time 
profile  in  the  limestone. 

In  addition  to  the  general  success  of  recording  directly  in  a  rock 
material,  there  are  several  other  noteworthy  aspects  of  this  particular 
test  record: 

(a)  Records  obtained  from  faster  sweeping  oscilloscopes  show  a  two 
step  rise,  the  first  corresponding  to  a  pressure  of  ~29  kbar  and  the 
second  to  a  peak  pressure  of  ''-88  kbar.  The  separation  between  waves  is 
<0. 1  /zsec. 


*  Sf«r|M  liasstoat,  p„  =  2.40  gfcc,  60*  fossil  frsfssats,  20*  biocsstic  dobris,  20*  cslcits;  aaslysis 
frow  hi.  20. 
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The  first  wave  had  not  been  detected  in  previous  work  with  Spergen 
limestone, 20  although  its  existence  was  expected.  Direct  comparison  of 
peak  pressure  results  with  previous  work  is  difficult  since  the  driver 
systems  were  not  the  same.  However,  the  peak  pressure  of  the  second  wave 
is  about  what  would  be  expected  from  the  limestone  Hugoniot  data  of  Gregson 
et  al. ,  and  the  cross  curve  of  the  explosive  driver,  Baratol.  High  sen¬ 
sitivity  oscilloscopes  (~1. 5  kbar  resolution)  showed  no  elastic  precursor 
wave,  which  is  in  agreement  with  the  data  obtained  by  optical  techniques 
on  similar  thickness  samples. 


1 - 

CURRENT  TURN  ON, 

CONDUCTIVITY  PROBE 

J _ 

1 - 

CURRENT  TURN  ON, 

MANGANIN  WIRE 

J _ 

TIME  —  Z /imc  /  major  division 

FIG.  14  RESPONSE  OF  MANGANIN  IN  LIMESTONE,  PEAK  PRESSURE  ~88  KBAR 

(b)  The  total  recording  duration  (before  wire  breakage)  was 
~17  /xsec.  Stretching  of  the  wire,  manifested  in  the  record  as  the  first 
increase  in  voltage  after  initial  peak,  occurred  after  6  /xsec  or  approxi¬ 
mately  at  the  arrival  time  at  the  wire  of  the  relief  waves  originating  at 
the  peripheral  gage-air  boundary.  In  effect,  the  wire  assembly  and  leads 
moved  in  pressure  and  particle  velocity  equilibrium  with  the  limestone, 
indicating  that  if  this  gage  system  were  in  an  all  limestone  medium,  the 
recording  duration  would  be  limited  by  the  length  of  the  limestone  gage 
assembly. 

A  similar  test  on  an  all 'tuff*  gage  gave  a  record  as  shown  in  Fig.  15. 
The  general  construction  was  similar  to  that  of  Fig.  13  and  is  shown  in 

•  NTS,  mi  v  tu,  p0  =  1.4  «/*c. 
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FIG.  15  RESPONSE  OF  MANGANIN  IN  TUFF,  PEAK  PRESSURE  ~  56  K BAR 


Fig.  16.  A  peak  pressure  of  56.4  ±  5  kbar  was  measured.  The  exact 
Hugoniot  for  this  sample  of  tuff  is  not  known;  however,  if  Lombard’s21 
curves  for  U  12 B  tuff  are  averaged,  a  pressure  of  60  1  6  kbar  would  be 
expected. 

The  records  from  two  granite*  gages  are  shown  in  Figs.  17  and  18. 

The  general  increase  in  noise  over  the  limestone  and  tuff  records  is  due 
to  the  granite  and  is  thought  to  be  a  piezoelectric  signal  from  the  in¬ 
cluded  quartz.  From  the  records  it  can  be  seen  that  at  the  higher  pressure, 
~280  kbar,  the  conductivity  of  granite  (measured  by  conductivity  probe) 
increases  (manifested  as  an  apparent  negative  pressure,  ~6  jitsec  after  peak) 
to  the  same  order  of  magnitude  as  that  of  the  Manganin  resistance  wire  and 
is  therefore  not  useful  as  a  gage  insulator.  These  results  have  prompted 
work  on  a  more  universal  gage  system  which  is  outlined  in  the  next  section. 

Although  the  peak  pressures  measured  by  the  wire  in  granite  were  ap¬ 
proximately  those  expected  from  the  granite  Hugoniot, 21  peak  pressures  were 
not  known  sufficiently  accurately  to  compare  the  Manganin  response  with 
Manganin  in  C-7  calibration. 


•  Herd  Hat  granite,  Are*  15  NTS,  pg  =  2,7  g/cc. 
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FIG.  17  RESPONSE  OF  MANGANIN  IN 


-ir 


FIG.  18  RESPONSE  OF  MANGANIN  IN 


D.  ADDITIONAL  GAGE  DEVELOPMENT 


The  ability 
rather  than  at  a  free  surface, 
to  field  tests.  The  preceding 


FIG.  19  MYLAR-MANGANIN-MYLAR 
FOIL  ASSEMBLY 


advantageous  to  laboratory  as  well  as 

with  Manganin  in  various  materials 
have  shown  that,  if  the  test  speci¬ 
men  is  an  insulator,  transducers 
can  be  constructed  from  them  and 
an  apparently  reliable  profile 
measured.  However,  should  the 
material  become  a  conductor  under 
shock  loading  the  transducer  ob¬ 
viously  fails.  Several  methods 
have  been  tried  to  insulate  the 
wire  from  its  surroundings  and  yet 
maintain  the  block  construction. 

In  all  trials,  aluminum  has  been 
used  as  the  test  material  because 
it  is  a  conductor  of  well-known 
equation  of  state,  thereby  permit¬ 
ting  simultaneous  measurements  of 
the  configuration  and  of  the 
Manganin  peak-pressure  calibration 
in  a  medium  other  than  C-7  epoxy. 


to  measure  pressure- time  profiles  within  a  material, 
i  s 

results 


The  configuration  shown  in  Fig.  19  was  found  to  be  satisfactory  for 
short  duration  recording,  <1  /isec.  It  is  constructed  by  compressing  the 
standard  four  lead,  0.003-inch-diameter  Manganin  wire  assembly  (see  Fig.  1) 
into  a  thin  ribbon,  0.001  inch  thick.  The  ribbon  is  placed  between 
~0. 0013-inch-thick  Mylar  films  and  the  assembly  compressed  between  heated 
(110°C)  anvils.  The  total  thickness  of  the  Mylar-Manganin-Mylar  assembly 
is  slightly  more  than  the  thickness  of  the  original  Manganin  wire  and,  ac¬ 
cordingly,  has  a  comparable  rise  time.  Although  Mylar  has  been  interposed 
in  the  system,  it  is  assumed  that  the  equilibrium  time  of  the  assembly  is 
sufficiently  short,  <0.1  jusec,  t0  cause  no  perturbations  in  the  shock 
profile. 

The  results  of  one  aluminum  block-foil  assembly  are  shown  in  Fig.  20. 
The  nature  of  the  failure  after  Ml. 8  nsec  has  not  yet  been  determined. 
However,  variations  of  this  basic  assembly  have  yielded  recording  durations 


TIME 


I  fisec/cm 


FIG.  20  RESPONSE  OF  MANGANIN-MYLAR  FOIL  IN  ALUMINUM, 
PEAK  PRESSURE  ~117  KBAR 


up  to  ~7  Msec,  about  the  limit  imposed  by  the  particular  overall  dimensions 
used.  The  fidelity  of  recording  after  peak  on  these  transducers  haB  not 
been  measured. 

E.  RECORDING  CIRCUITS 

Standard  high  frequency  pulse  techniques  have  been  observed  in  the 
recording  circuits  shown  in  Fig.  4.  The  change  in  Manganin  resistance  is 
measured  by  a  four- terminal  system,  two  voltage  leads  and  two  current  leads. 
The  use  of  a  four-lead  system  minimizes  errors  introduced  by  the  possible 
change  in  resistance  at  the  junction  of  the  wire  and  current  posts.  Ef¬ 
fectively,  the  measured  change  in  resistance  is  that  between  the  voltage 
taps,  which  permits  measurements  to  be  confined  to  the  center  region  of 
the  wire  only;  an  example  is  the  case  of  the  reflected  wave  shots  of 
Section  IV  B. 

A  constant  current  power  supply  (Fig.  21)  minimizes  errors  due  to 
stray  inductance.  The  supply  is  triggered  just  prior  to  shock  arrival  at 
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the  wire  to  allow  display  of  the  current  turn-on  signal.  This  signal  ap¬ 
pears  as  a  voltage  step,  Vfl,  which  is  proportional  to  the  initial  wire 


resistance,  fi0,  and  is  followed  by  a  voltage  step,  AF,  proportional  to 


the  change  in  resistance  Aft,  at  the  time  of  shock  arrival.  With  this 
method,  the~  necessity  of  voltage  and  current  calibration  is  eliminated. 
The  actual  quantity  measured,  AV,/Kfl,  is  related  to  Aft/ft.  by 


A ft  AF  „ 

=  —  +  /3 

*o 


(13) 


where  /3  accounts  for  the  slight  change  in  wire  current  due  to  the  parallel 
resistances  of  the  differential  voltage  circuits.  Values  of  /?  vs.  AF/F# 
are  shown  in  Fig.  22.  For  field  use  the  correction  is  insignificant  but 
must  be  considered  in  the  peak  pressure  calibration  discussed  in 
Section  IV  B. 


FIG.  22  RECORDING  CORRECTION  FACTOR  /S,  vs.  AV/V0 


The  differential  recording  system  provides  the  normal  benefits  of 
common  mode  rejection,  thereby  assuring  measurement  of  only  the  change 
in  resistance  of  the  center  portion  of  the  wire.  However,  care  must  be 
taken  in  the  system  balance;  for  example,  in  the  voltage  dividers  formed 
by  the  voltage  tap  series  and  termination  resistors.  Although  it  is  not  a 
problem  in  laboratory  tests,  field  tests  requi re  more  careful  cable  balance. 


F.  FIELD  TESTS 


1.  General 

Early  field  tests  are  reported  in  Reference  4.  More  recent  tests 
have  been  conducted  in  conjunction  with  the  FLAT  TOP  series  of  events  to 
provide  fielding  experience  with  various  gage  and  trigger  systems.  All 
transducers  were  of  the  Manganin  in  C-7  configuration.  Since  the  Hugoniot 
of  the  surrounding  medium  was  generally  not  known,  the  recorded  pressure 
profiles  have  not  been  related  to  the  surrounding  medium.  The  transducer 
results  of  each  event  are  described  below  in  chronological  order.  Re¬ 
cording  and  overall  performance  results  are  described  in  another  report. 2 

2.  Flat  Top  II 

In  this  test  conducted  in  playa,  a  50-11  transducer4  using  a  bridge 
circuit4  was  placed  ~6  inches  from  the  source  (TNT)-sand  interface.  Since 
the  shock  tilt  and  shock  properties  of  the  reconstructed  medium  in  the  im¬ 
mediate  vicinity  of  the  gage  were  not  known,  the  pressure  amplitude  results 
are  of  dubious  value.  However,  the  function  tests  of  the  gage  and  recording 
systems  were  successful. 

The  records  are  shown  in  Fig.  23.  Curve  “  a”  is  a  reproduction  of  the 
actual  trace  recorded,  in  which  the  pressure  pulse  occurred  just  at  the  end 
of  the  power  supply  turn-on  signal.  Because  the  gage  signal  appeared  early, 
it  can  be  concluded  that  there  was  considerable  shock  tilt  and  that  the 
shock  wave  reached  the  PZT  piezoelectric  trigger  probe,  placed  three  inches 
in  front  and  on  top  of  the  gage,  just  shortly  prior  to  reaching  the  gage. 
These  conclusions  are  substantiated  by  the  relatively  slow  rise  time  of  the 
gage  signal.  The  "true”  recorded  pressure- time  profile,  shown  as  curve 
“  c,”  can  be  obtained  from  the  baseline  shift  during  this  interval,  curve 
“b”;  fortunately,  this  was  recorded  prior  to  the  shot.  The  peak  pressure 
recorded  was  30.4  ±  1  kbar.  Approximately  7.5  jzsec  after  shock  arrival, 
the  gage  wire  began  to  stretch  due  to  the  peripherally  originating  relief 
waves.  After  ~30  fJ.se c,  the  wire  was  seen  to  break.  The  mismatch  between 
gage  and  soil  is  ntrt  known,  as  the  soil  was  hand-tamped  around  the  gage. 

It  is  noteworthy  that  the  gage  remained  intact  in  this  environment  for  a 
considerably  longer  time  than  a  similar  shot  in  air;  i.e.,  although  the 
arrival  time  of  the  relief  wave  is  nearly  the  same,  the  amplitude,  as 
might  be  expected,  is  considerably  less.  As  pointed  out  previously,  the 
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gage  is  measuring  a  profile  which  exists  in  the  gage  insulator  and  which 
departs  from  the  profile  in  the  surrounding  medium  at  some  time  after 
peak  pressure,  depending  upon  the  gage  dimensions. 

3.  Flat  Top  III 

a.  GAGE  CONFIGURATION  AND  PLACEMENT 

Four  Manganin  wire  gages,  three  1-H  differential  and  one  50-fl, 
were  placed  abutting  native  playa  at  distances  of  53,  60.5,  73,  and 
94.5  inches  from  shot  center.  Closest  distances  to  a  TNT  surface  turned 
out  to  be  on  the  order  of  2,  6,  26,  and  42  inches,  respectively.  The 
general  placement  configuration  is  shown  in  Fig.  24(a)  and  (b),  the  re¬ 
cording  schematic  in  Fig.  25,  and  a  typical  1-fl  differential  gage  in 
Fig.  26.  This  gage  was  constructed  to  provide  a  20-fisec  duration  re¬ 
cording  at  100  kbar  and  was  identical,  except  in  diameter  and  thickness, 
to  standard  laboratory  gages.  The  50-0  gage  was  a  sel f- triggering  unit 
constructed  for  a  previously  scheduled  NTS  test.  The  power  supply  cir¬ 
cuitry  for  each  type  has  been  described  in  detail  elsewhere. 4 

.6.  RESULTS 

Records  obtained  from  gage  numbers  EP  1  and  EP  4  (see  Fig.  25), 
both  standard  1-0  gages,  are  shown  in  Figs.  27  and  28.  In  both  cases  the 
pressure  signal  occurred  during  power  supply  turn  on;  the  waveform  shown 
is  a  difference  signal  between  the  test  trace  recorded  prior  to  the  event 
and  the  event  signal.  The  gages  were  positioned  prior  to  HE  stacking  and 
were  referred  to  shot  center;  i.e.,  the  trigger  probes,  PZT  piezoelectric 
ceramic  elements,  were  placed  ~3  inches  in  front  of  the  gages  on  a  sphere 
radius  as  shown  in  Fig.  24(a).  Therefore,  the  early  arrival  of  the  shock 
wave  at  EP  1,  with  reference  to  its  trigger  probe,  is  not  surprising. 
Similarly,  it  can  be  expected  that  the  pressure  profile  would  be  rather 
unlike  a  plane  wave  in  the  vicinity  of  EP  1,  and  the  shocks  traveling  per¬ 
pendicular  to  the  gage  axis  would  distort  the  signal.  The  profile  recorded 
by  EP  4  indicates  an  initially  lower  pressure  than  expected  (25  1  10  kbar 
vs.  50  kbar)  but  a  final  pressure,  before  wire  breaking,  of  56  ±  10  kbar. 
The  step  increases  may  be  due  to  multiple  shocks  from  the  various  layers 
of  HE.  Shock  arrival  time,  with  reference  to  gage  trigger  probe,  is  very 
close  to  that  expected,  and  the  intermixing  of  event  signal  with  power 
supply  is  caused  by  the  unusually  long  turn-on  signal  associated  with  this 
single  cable  differential  system. 
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FIG.  28  1-0  DIFFERENTIAL  GAGE  RECORD,  EP4,  FLAT  TOP  III 


The  50-0  gage  (EP  2)  records  are  off-scale  on  all  oscilloscopes. 
It  is  not  known  whether  this  is  trigger  noise  coupled  from  failure  of  the 
silicon  controlled  rectifier  in  the  gage  or  an  unusually  large  signal.  In 
view  of  the  pressure  amplitudes  recorded  by  EP  1  and  EP  4,  the  former 
possibility  is  more  likely.  No  oscilloscope  traces  were  recorded  from  the 
last  gage,  EP  5.  The  reason  for  this  is  not  known,  but  it  is  suspected 
that  the  delay  oscilloscope  (see  Fig.  25)  must  have  triggered  prematurely 
since  gage  system  EP  4,  from  which  the  trigger  was  derived,  functioned 
satisfactorily. 

4.  Flat  Top  I 

Four  gages  of  the  type  shown  in  Fig.  26  were  constructed  for  this 
event.  Two  were  mounted  on  8-inch  diameter  by  ~18-inch-long  marble  cy¬ 
linders  which  were  placed  in  the  limestone  shot  medium.  One  was  placed 
in  a  column  of  grout  and  one  on  the  HE  stack.  (Details  are  described  in 
Reference  22.)  Records  were  obtained  from  the  latter  two  gage  systems. 

Of  the  other  two  gages,  one  functioned  satisfactorily  except  that  the 
noise  level  at  the  time  of  recording  masked  the  pressure  record  and  one 
recording  system  appeared  to  have  triggered  early. 

The  record  from  the  gage  on  the  HE  stack  (reduced  as  in  the  pre¬ 
ceding  shots)  is  shown  in  Fig.  29.  Both  the  duration,  ^10  jusec,  and  the 
peak  pressure,  ~95  kbar,  are  lower  than  expected.  One  postshot  test  was 
conducted  in  which  a  gage  of  identical  construction  was  subjected  to  a 
shock  from  a  TNT  explosive  pad.  The  predicted  peak  pressure  was  145  to 
150  kbar,  and  a  value  of  147  kbar  was  recorded,  the  duration  before 
stretching  being  as  expected,  ~20  /xsec.  Since  the  only  major  difference 
was  that  the  laboratory  gage  was  cemented  to  the  TNT  whereas  the  field 
gage  was  not,  the  coupling  of  gage  to  the  TNT  is  suspect. 
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FIG.  29  1-ft  DIFFERENTIAL  GAGE  RECORD,  HE  STACK,  FLAT  TOP  I 

The  grout  gage  record,  again  reduced,  is  shown  in  Fig.  30.  The  peak 
pressure  reached  is  ~34  kbar  and  the  duration  at  peak  pressure  is  ~20  /zsec 
No  further  analysis  has  been  made  since  the  relation  between  the  profile 
in  grout  and  that  in  the  surrounding  limestone  is  not  known. 


1-0  DIFFERENTIAL  GAGE  RECORD,  GROUT  COLUMN,  FLAT  TOP  I 


V  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

The  general  fidelity  and  response  of  a  Manganin  in  C*7  epoxy  trans¬ 
ducer  system  appears  suitable  for  use  as  a  field  test  gage  i'n  a  pressure 
range  from  as  low  as  field  recording  limitations  will  permit,  ~ 5  kbar, 
to  at  least  200  kbar  gage  pressure.  Because  C-7  is  a  material  of  low 
shock  impedance  when  compared  to  most  rocks,  it  can  be  used  to  measure 
pressure  profiles  of  considerably  higher  rock  pressures.  The  pressure 
limit  in  soils  is  approximately  that  of  the  gage  since  C-7  shock  imped¬ 
ance  is  very  similar  to  tuff,  .alluvium,  and  playa.  In  rocks  where  the 
gage  and  medium  are  mismatched,'  the  pressure  profile  existing  in  the 
medium  can  be  obtained  if  its  equation  of  state  is  known.  In  the  corre¬ 
sponding  laboratory  example  chosen  for  analysis,  C-7  epoxy  on  aluminum, 
the  profile  can  be  obtained,  accurate  to  within  10  percent,  by  a  point 
for  point  impedance  match  solution;  i.e.,  the  aluminum  pressure  can  be 
obtained  from  the  appropriate  Hugoniots  and  cross  curves. 

The  recording  duration  of  the  present  type  of  gage,  when  used  in 
media  of  dissimilar  shock  characteristics,  is  necessarily  limited  by 
gage  dimensions.  At  moderate  pressures,  <  100  kbar,  and  slowly  decaying 
waves,  laboratory  tests  have  shown  this  to  be  ^20  fi sec  for  a  6-inch 
diameter  by  4-inch-thick  gage.  Long  duration  waves  at  higher  pressures 
in  the  low  impedance  C-7  are  beyond  the  test  capability  established  by 
current  laboratory  explosive  size  limitations.  Recording  directly  in 
higher  impedance  media  may  remove  this  test  restriction.  However,  higher 
pressures  will  require  either  that  there  be  a  better  match  between  the 
gage  insulators  and  gage  cur rent- vol tage  leads  or  that  the  shock  response 
of  thin  foils  leads  be  better  understood.  This  latter  problem  is  a  de¬ 
parture  from  the  one-dimensional  systems  considered  in  the  present  work. 
However,  such  a  two-dimensional  study  has  important  application  to  field 
test  gage  placement,  e.g.,  the  behavior  of  a  column  of  grout  in  a  test 
medium  of  differing  Hugoniot  characteristics.  If  the  profile  in  the 
medium  can  be  obtained  from  measurements  in  the  grout,  the  general  gage 
placement  is  greatly  simplified. 
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The  feasibility  of  recording  directly  in  other  media  has  been  demon* 
strated.  However,  the  general  response  of  Manganin  in  these  media  must 
be  determined.  In  addition,  a  more  universal  wire  configuration  is  re* 
quired  to  permit  recording  in  materials  other  than  insulators.' 

Field  test  results  are  difficult  to  summarize  since  gage  placement 
geometry  differed  considerably  from  the  one-dimensional  systems  for  which 
the  gage  wa$  laboratory  tested.  Hence,  little  can  be  concluded  from  peak 
pressure  measurements  about  gage  functioning.  Durations  of  recording, 
in  most  cases,  were  about  those  expected. 

B.  RECOMMENDATIONS 

1.  Manganin  in  C-7  Epoxy 

The  Manganin  in  C-7  gage  has  been  tested  over  part  of  the  range  of 
desired  gage  response  characteristics  and  been  found  to  perform  satis¬ 
factorily.  Since  long  duration  recording  at  high  pressures  is  a  goal, 
it  appears  that  response  under  these  conditions  muat  be  either  inferred 
from  lower  pressure  work  or  measured  in  actual  field  tests.  Although 
the  overall  response  is  satisfactory  and  is  approaching  a  working  tool 
for  certain  ranges  of  pressures  and  recording  times,  it  is  recommended 
that  the  limitations  of  the  gage  be  determined  by: 

a)  completing  calculations  of  wire  temperatures  due 
to  conduction  of  heat  from  the  epoxy, 

b)  completing  previously  proposed  work  on  lead- 
insulator  matching  and  subsequent  testing  at 
high  pressures,  and 

c)  measuring  hysteresis  and  duration  of  recording 
in  a  similar  pressure  range. 

2.  New  Gages 

The  most  promising  recent  development  is  the  direct  measurement 
technique  outlined  in  Section  IV  D.  The  eventual  success  of  this  method 
depends  on  continuing  several  areas  of  work: 

a)  developing  methods  of  measuring  pressure-time 
profiles  with  Manganin  foils  in  conductors, 

b)  measuring  Manganin  response  in  various  media 
for  comparison  with  that  obtained  in  C-7  and, 

in  particular,  measuring  response  to  anisotropic 
stress. 
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3.  Two-Dimensional  Shock  Studies 


Concurrent  with  the  above  work,  it  is  recommended  that  a  program  be 
initiated  to  investigate  the  behavior  of  materials  under  oblique,  nonsteady 
state  shocks  by  examining  the  relation  between  stress  profiles  in  adjacent 
media  of  different  equations  of  state.  Although  some  theoretical  work  is 
possible,  the  emphasis  would  be  on  experimental  work.  The  goals  of  this 
study  should  be  1)  to  predict  gage  component  (lead)  requirements  and 
2)  to  provide  data  useful  in  field  installation  of  gages. 
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APPENDIX  A 


DETERMINATION  OF  TEST  MATERIAL  PRESSURE 
PROFILE  FROM  TRANSDUCER  RECORD 
BY  METHOD  OF  CHARACTERISTICS 

Consider  the  problem  of  an  aluminum  plate  in  contact  with  a  slab  of 
C-7  with  a  resistance  wire  gage  imbedded  in  the  C-7  two  mm  from  the  inter¬ 
face.  At  the  back  of  the  aluminum  plate  a  shock  is  induced  which  passes 
through  the  aluminum  and  into  the  C-7.  As  the  shock  passes  the  gage,  a 
pressure-time  history  is  recorded.  Reading  values  of  pressure  at  arbitrary 
time  intervals,  0.4  /xsec,  produces  a  series  of  pressure-time  steps  which 
can  be  used  to  generate  the  solution. 

A  working  model  of  the  transformation  in  the  distance-time  (*, t) 
plane  is  required.  This  depicts. qualitatively  what  happens  during  the 
shock  process,  Fig.  A-l  serving  as  the  model. 

The  construction  of  the  model  is  based  on  the  simple  qualitative 
attributes  of  shock  waves,  such  as  the  production  of  both'  a  reflected 
wave  and  a  transmitted  wave  at  the  interface  between  materials  of  dif¬ 
ferent  impedances.  From  this,  the  model  is  constructed  as  follows: 


FIG.  A-l  REFLECTION  OF  A  DECAYING  SHOCK  AT  ALUMINUM-EPOXY  INTERFACE 
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The  initial  shock  in  the  aluminum  is  drawn  to  the  interface  between  the 
aluminum  and  the  C-7.  (The  slope  is  arbitrary  at  this  time.  When  the 
solution  is  actually  being  generated,  it  will  depend  on  the  pressure  in 
the  aluminum. )  Since  there  is  an  impedance  mismatch,  there  will  be  both 
reflected  and  transmitted  waves.  By  using  the  method  of  impedance  matching, 
it  can  be  seen  in  Fig.  A-2  that  a  shock  is  transmitted  to  the  C-7  (line  3 
in  Fig.  A- 1 )  and  a  rarefaction  is  reflected  into  the  aluminum.  This  rare¬ 
faction  is  depicted  in  Fig.  A-l  in  its  usual  form,  i.e.,  a  rarefaction 
"fan”  (lines  2). 


FIG.  A-2  HUGONIOT  AND  CROSS  CURVES  IN  THE  PRESSURE-PARTICLE 
VELOCITY  PLANE 


When  the  shock  initially  stikes  the  back  (see  Fig.  A-l)  of  the 
aluminum  plate,  this  region  of  the  plate  accelerates  to  some  nonzero 
particle  velocity.  Similarly,  once  the  shock  reaches  the  interface  and 
the  resistance  wire,  both  assume  some  particle  velocity.  The  wire  veloc¬ 
ity  is  assumed  to  be  that  of  the  moving  mass  of  C-7,  that  is,  the  same 
as  the  interface  velocity.  The  interface  line  (4),  the  line  representing 
the  back  of  the  plate  (5),  and  the  relief  waves  originating  at  the  back 
of  the  aluminum  plate  are  shown  schematically  in  Fig.  A-l.  In  accordance 
with  the  theory  of  propagation  of  small  pressure  pulses,  the  slope  of 
each  of  these  right-traveling  waves  will  be  dx/dt  =  (u  +  c)  where  u  and 
c  are  the  average  particle  velocity  and  sound  velocity  in  the  two 
regions  separated  by  the  relief  wave.  If  the  relief  waves  were 
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left-traveling,  their  slope  would  be  dx/dt  =  u  -  c.  Note  that  when  each 
relief  wave  reaches  the  interface,  a  compression  is  reflected  into  the 
aluminum  and  a  relief  wave  is  transmitted  into  the  C-7.  The  number  of 
relief  waves  drawn  in  the  final  solution  is  dependent  on  the  number  of 
pressure- time  points  that  are  read  from  the  gage  record. 

The  numerical  solution  can  now  be  Started.  First,  the  portion  of 
the  model  around  the  shock  is  enlarged,  and  each  region  (that  is,  each 
bounded  area)  is  labeled  (Fig.  A-3). 


As  an  approximation,  it  is  assumed  that  in  each  region  the  particle 
velocity  and  the  sound  velocity  are  constant  throughout  the  region— the 
smaller  the  region,  the  more  nearly  correct  this  approximation.  With 
this  assumption,  each  region  of  the  x,t  plane  will  map  into  a  single 
point  on  a  plot  of  particle  velocity  vs.  sound  velocity.  Thus  the 
introduction  of  the  particle  velocity-sound  velocity  (u,c)  plane  becomes 
advantageous. 

Without  discussing  the  theory  of  characteristics,  it  can  be  shown1 
that  for  the  Murnaghan  equation  of  state  (entropy  assumed  constant) 
p  =  A [{p/p9 )*-l] ,  the  equations  of  the  two  families  of  curves  that 


G.  It  Fovles,  "Attaaaatioa  of  a  Shock  Vara  Prodacad  ia  a  Solid  by  a  Flyiaf  Plata  *  J.  Appl.  ftya.  31, 
6S$>  661  (I960) 
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represent  characteristics  in  the  u, c  plane  are 


„  *  dc  (if  -  1) 

C  +  characteristics:  -  =  - 

du  2 

„  .  .  .  dc  (K  -  1) 

C  ~  characteristics:  -  =  -  - — 

du  2 

where  K  is  the  constant  in  the  given  Murnaghan  equation  Thus  the 
characteristics  are  straight  lines  in  the  u,c  plane. 

In  any  single  material,  regions  in  the  x,t  plane  which  are  separated 
only  by  left-traveling  waves  become  points  in  the  u,c  plane  which  lie 
along  common  C~  characteristics;  for  instance,  in  Fig.  A-3,  regions  a, 
b,  c,  and  d  lie  along  a  common  C~  characteristic  (Fig.  A-4). 


(A-l) 

( A-  2 ) 


FIG.  A-4  CHARACTERISTICS  IN  THE  (NORMALIZED)  SOUND  VELOCITY -PARTICLE 
VELOCITY  PLANE.  It  is  common  to  plot  c/cq  vs.  u/c0  whors  cQ  is  tho  sound 
velocity  in  a  reference  region  labeled  "0”  in  Figure  4. 

Similarly,  regions  separated  by  only  right-traveling  waves  lie  along 
common  C +  characteristics.  Using  the  approximation  that  the  shock  is  an 
isentropic  compression,  regions  0  and  a  in  Fig.  A-3  are  shown  in  Fig.  A-4 
lying  on  a  common  C *  characteristic.  Note  that  regions  0  and  b,  sepa¬ 
rated  by  both  left-  and  right-traveling  waves,  do  not  lie  on  any  common 
characteristic.  It  is  clear  that  if  u  and  c  are  known  in  region  d,  for 
example,  and  either  u  or  c  is  known  in  regions  a,  6,  and  c,  then  all  four 
points  (a,  b,  c,  d)  can  be  located  in  the  u,c  plane,  since  they  must  lie 
along  a  common  C~  characteristic.  Once  both  sound  velocity  and  particle 
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velocity  are  known  at  each  of  these  points,  the  lines  which  represent 
the  waves  separating  regions  a  and  b.,  b  and  c,  and  c  and  d  can  be  drawn. 

This  is  possible  since,  as  mentioned  previously,  small  pressure  pulses 
propogate  with  Velocities  dx/dt  =  u  +  c  for  right-traveling,  u  -  c  for 
left- traveling  waves. 

Values  of  sound  velocity,  c,  can  be  obtained  as  follows:  There 
exists  a  one-to-one  correspondence  between  sound  velocity  [defined  by 
c2  =  dp / dp]  and  pressure  under  the  assumption  of  constant  entropy.  In 
the  case  of  a  Murnaghan  equation  of  state, 

p  =  A[(p/p0)K  -  1]  (A- 3) 

and 

c2  =  UK/p0)(p/p0)l‘-1  =  K(p  +  A)/p  .  (A-4) 

To  clarify  the  above  general  discussion,  the  actual  calculations 
for  a  specific  problem  will  be  partially  reproduced.  The  pressure  gage 
readings  are  shown  in  Fig.  A-5.  From  the  initial  reading,  the  pressure 
in  region  e  (Fig.  A- 3 )  of  the  C-7  is  148  kbar;  using  the  Murnaghan  equa¬ 
tion  of  state  for  aluminum,  this  implies  a  particle  velocity  of  2.16  mm//isec 
and  a  shock  velocity  of  6.09  mm/psec  in  the  C-7.  Since  particle  velocity 
and  pressure  must  be  continuous  across  the  interface,  in  region  d  of  the 
aluminum,  p  3  148  kbar  and  u  3  2.16  mm/psec.  After  the  shock  in  the 
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C«7  is  drawn  to  proper  slope  ( dx/dt  =  6.09  mm/^tsec)  and  the  interface  to 
its  slope  ( dx/dt  =  =  2.16  mm//i,sec),  d  is  located  in  the  u,c  plane  and 

the  C-  characteristic  [slope  =  -(K  ~  l)/2]  is  drawn  through  it  [Fig.  A-6(a)] 
Points  a,  b,  and  c  in  the  u,c  plane  must  lie  on  the  characteristic  through 
d.  The  location  of  point  a  in  the  u,c  plane  can  be  found  by  either  (1)  the 


(NOT  TO  SCALE) 
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FIG.  A-6 

CHARACTERISTICS  CURVES 

impedance  match  method  or  (2)  drawing  the  C+  characteristic  through 
point  0  [slope  ( K  ~  l)/2]  to  intersect  the  C~  characteristic  through 
d.  The  intersection  must  be  point  a,  as  was  shown  in  the  discussion 
of  the  model. 

By  using  the  impedance  match  procedure,  the  pressure  in  region  a 
is  303  kbar.  This  corresponds  to  U f  =  7.31  mm/fj.sec  and  up  -  1 . 48  mm/fisec . 
Since  region  a  is  assumed  to  extend  to  the  back  of  the  aluminum  plate 
(to  be  discussed  later),  the  velocity  of  the  back  of  the  aluminum  plate 
is  just  up  in  region  a,  or  dx/dt  =  1.48  mm/pLsec.  Between  regions  a  and 
d  there  is  the  rarefaction  fan;  the  pressure  drop  across  each  of  the 
elements  of  the  fan  is  assumed  equal,  hence  the  pressures  are  p ^  =  303  kbar, 
p b  =  252,  pc  =  200,  pd  =  148.  Since  sound  speed  is  a  unique  function  of 
pressure,  this  locates  points  fa  and  c  in  the  u,  c  plane.  Thus  the  waves 
separating  regions  a,  fa,  c,  and  d  can  De  drawn  with  dx/dt  =  u  -  c  in  each 
case.  In  Fig.  A-6(b), 

For  wave  (1),  dx/dt  =  u  -  c  =  -5.89 

For  wave  (2),  dx/dt  =  u  -  c  =  -5.34 

For  wave  (3),  dx/dt  =  u  -  c  -4.78 


This  calculation  completes  the  first*' row*'  of  regions  in  the  aluminum. 
To  begin  the  next  row,  containing  regions  g  through  k  in  Fig.  A-6(b),  the 
first  pressure  drop  measured  from  the  gage  is  determined.  In  this  case 
the  drop  occurs  0.4  fiae c  after  the  initial  shock  reaches  the  wire.  In 
Fig.  A-6(b)  the  pressure  in  region  /  in  the  C-7  is  140  kbar  (from  gage 
record),  which  corresponds  to  a  particle  velocity  of  2.08  mm//zsec.  The 
wave  is  drawn  with  this  slope  through  the  point  on  the  gage  line  0.4  /usee 
above  the  point  at  which  the  initial  shock  reached  the  wire.  Again, 
pressure  and  particle  velocity  must  be  continuous  across  the  interface, 
so  the  location  of  point  g  in  the  u,c  plane  is  known  [Fig.  A-6(a)]. 

Drawing  the  C-  characteristic  through  g  in  the  u,c  plane  yields  the 
line  on  which  points  h,  i,  j,  and  k  must  lie.  But  d  and  h  are  separated 
only  by  a  right- travel ing  wave;  hence  they  lie  in  a  common  C+  character¬ 
istic  in  the  u,c  plane.  The  intersection  of  the  C+  characteristic  through 
d  and  the  C-  characteristic  through  g  is  h.  Similarly,  points  i,  j,  and 
k  are  obtained  by  drawing  the  C+  characteristics  through  c,  b,  and  a, 
respectively,  until  they  intersect  the  C-  characteristic  through  g. 
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The  lower  particle  velocity  in  /  (as  compared  to  e)  indicates  that 
the  slopes  of  the  interface  and  wire  lines  change;  likewise  the  lower 
particle  velocity  in  k  (as  compared  too)  implies  that  the  slope  of  the 
back  of  the  aluminum  plate  changes. 

The  next  “row"  of  regions  is  calculated  in  exactly  the  same  manner. 

The  state  in  C-7  is  known  from  the  gage  record;  this  state  defines  the 
slope  of  the  relief  wave.  The  state  in  the  aluminum  at  the  interface  is 
required  to  have  the  same  pressure  and  particle  velocity  as  the  state  in 
C-7  though  the  sound  velocity  is  different  This  then  locates  the  specific 
C-  characteristic  on  which  all  points  in  the  row  lie;  by  extending  the  C+ 
characteristics  through  the  points  d  and  h,  c  and  i,  b  and  j,  etc.,  all 
the  points  in  the  row  are  determined.  The  slopes  of  the  waves  in  the 
x,  t  plane  are  then  calculated  as  specified  above.  This  procedure  is  con¬ 
tinued  until  the  data  from  the  pressure  gage  record  are  exhausted.  The 
values  of  the  flow  variables  in  each  region  are  given  in  Table  A-l. 


Table  A-l 
FLO*  PARAMETERS 


iUfioa 

Praaaara 

(kbar) 

(aa/wte ) 

II 

(n^taae) 

a 

303 

7.65 

mm 

b 

252 

7.34 

c 

7.01 

BK 

d 

148 

6.65 

BBS 

e 

148 

6.38 

B9 

f 

140 

6.26 

2.08 

t 

140 

6.59 

2.08 

h 

136 

6.56 

2.10 

i 

*"  187 

6.92 

1.87 

j 

236 

7.24 

1.66 

h 

292 

7.58 

1.45 

Reconstructing  the  undistorted  wave  in  aluminum  by  this  procedure 
necessarily  limits  the  analysis  to  pressures  equal  to  or  less  than  the 
peak  as  recorded  by  the  gage.  This  is  because  any  relief  waves  which 
have  intersected  the  initial  shock  before  it  reaches  the  gage  must  be 
neglected;  they  cannot  be  determined  from  the  data  available.  Thus  the 
pressure  in  region  a  of  Fig.  A-6,  which  was  found  to  be  303  kbar,  is 
almost  certainly  not  the  pressure  at  the  back  of  the  plate.  As  an  example 
of  this,  consider  the  relief  wave  drawn  as  a  dashed  line  in  Fig.  A-7;  this 
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intersects  the  shock  before  it  reaches  the  gage.  The  pressure  in  region 
z,  which  extends  to  the  back  of  the  aluminum  plate,  is  greater  than  that 
in  region  a.  But  since  there  is  no  way  by  which  the  existence  of  z  can 
be  known  from  data  at  the  gage,  it  must  be  assumed  that  region  a  extends 
to  the  back  of  the  plate.  Since  the  greater  pressure  in  region  z  corre¬ 
sponds  to  a  higher  shock  velocity  than  that  in  region  a,  the  true  profile 
of  the  shock  will  be  altered;  the  shock  including  the  effect  of  z  is 
labeled  “true  shock”  in  Fig.  A-7. 

The  results  of  the  analysis  are  given  in  Sec.  III-B. 


FIG.  A-7  CHARACTERISTICS  PATTERN  FOR  SHOCK  WAVE  DECAYING 
PRIOR  TO  REFLECTION  AT  ALUMINUM-EPOXY  INTERFACE 
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HUGONIOT  TEMPERATURE  CALCULATIONS  FOR  C-7  EPOXY 

By 

P.  H.  Moravek 
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APPENDIX  B 


HUGONIOT  TEMPERATURE  CALCULATIONS  FOR  C-7  EPOXY 


This  appendix  deals  with  calculations  of  shock  temperatures  reached 
in  the  C-7  insulator  due  to  shock  waves  of  amplitude  greater  than  500  kbar. 
The  Hugoniot  is  represented  by  the  equation1 


P„  =  78.16m  +  195.6m2  +  221.4 m3  kbar 


(B- 1) 


where 


G-HH 


(B-2) 


The  symbols  are  defined  below. 

A.  SYMBOLS  AND  CONSTANTS 

1.  Definitions  of  Symbols  Used  in  Text 

Cp  specific  heat  at  constant  pressure 

C^Q  initial- value  specific  heat  for  C-7,  cal/g°K 

<^,'0  4.186  x  10_2C^0,  initial-value  specific  heat,  kbar  cnfJV’  K 

Cy  specific  heat  at  constant  volume 

Cyj  specific  heat  in  region  I  for  300  K  <  T  <  902°K,  kbar  cm^g°K 

CytI  specific  heat  in  region  II  for  902°K  £  T,  kbar  cm^g°K 

Eg  specific  internal  energy  for  Hugoniot  states 

Eg  specific  internal  energy  for  isentropes 

I.(a)  (t  *  1,2,  3, 4)  (a  =  A,B )  integral  expression  depending  on 
each  Gruneisen  ratio  used 
P  pressure,  kbar 

Pg  Hugoniot  pressure,  kbar 

Pg  isentrope  pressure,  kbar 

S  entropy,  cal/g°K 

O 

T  temperature,  K 


1  John  0.  Erkaaa,  private  coMnnicatioa.  Since  Eq.  (B-l)  nan  fitted  to  data  only  ap  to  eboet  100  kbar, 
*e  ahoald  be  critical  of  reaalta  for  high  preaaoraa. 


73 


initial  temperature,  °K 
temperature  along  the  Hugoniot,  °K 
temperature  on  isentrope,  °K 
volume,  cm  3/g 

initial  vo]  e  at  P  =  0,  cm3/g 
'  ■  relative  volume  parameter 

a  constant  in  specific  heat  equation 
a  constant  in  specific  heat  equation 
Griineisen’s  ratio 

initial  value  of  Griineisen's  ratio  for  the  data 
available  on  C- 7 

initial  value  of  Griineisen's  ratio  for  Lucite 
(i  =  l,2,3,4)(a  =  A,B)  Griineisen’s  ratio  for  the 
cases  considered 


(V-‘)  ■ 


relative  volume  parameter 


upper  probable  error  in  measurement  of  T 
lower  probable  error  in  measurement  of  f 


2.  Values  of  Some  Constants 


-  0.387  cal/ g°K  at  T0„ 

«  1.855  x  10"  2  kbar  cm3/g°K 

»  3.63  x  10”2  kbar  cm3/g°K  -  6.00  cal/atom°K 
based  on  average  molecular  weight  326.2  g/mole 
and  average  atomic  weight  6.91  g/atom 

■  300°K 

-  0.8475  cm 3/ g 

■  0.106 
■  0.352 
=  0.790 
=  1.025 


As  reported  previously,  approximate  temperatures  of  the  C-7  insula* 
tor  were  determined  by  calculations  based  on  elastic  and  thermal  data  for 
Lucite  along  with  certain  simplifying  assumptions  such  as:  specific  heat, 
Cy ,  is  constant;  Griineisen’s  ratio, 


—  (dP\ 
Cy  \ST/ 


(B-3) 


P.  H.  Moravek,  Appendix  A,  SRI  Project  GPU* 37 13,  Progress  Report  No.  27,  December  17,  1963. 
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is  a  constant;  and  T/V  is  constant  for  a  second  calculation.  Lucite  data 
for  specific  heat,  compressibility,  thermal  expansion,  density,  and 
Hugonoit  pressures  were  used  because  some  of  the  properties  are  quite 
comparable  to  those  of  C-7,  or  at  least  were  believed  to  be  so;  also, 
considerably  more  data  were  available  for  Lucite.  There  are,  however, 
differences  between  C-7  and  Lucite  and  a  more  accurate  set  of  approxima¬ 
tions  is  given  in  the  present  report  to  account  for  the  C-7  Hugoniot  data, 
the  approximate  variation  of  C-7  specific  heat  with  temperature,  and  a 
range  of  straight-line  approximations  to  Griineisen's  ratio  as  a  function 
of  volume  based  on  experiment.  The  Griineisen  ratio  for  materials  is  really 
a  function  of  both  volume  and  temperature,  but  the  temperature  dependence 
is  believed  to  be  insigni ficant3 except  at  low  temperatures. 

The  previous  approximate  calculations  were  made  for  the  shock-pressure 
range  up  to  125  kbar  at  which,  using  either  T  constant  or  F/V  constant, 
the  temperature  becomes  of  the  order  of  10*  °K.  The  results  show  the  in¬ 
teresting  fact  that  at  125  kbar  the  temperatures  calculated  for  the  two 
cases  differ  by  less  than  seven  percent. 

B.  PROCEDURE  AND  RESULTS 

This  appendix  gives  calculated  Hugoniot  temperatures  for  pressures 
exceeding  500  kbar  for  four  slightly  different  Griineisen’s  ratios. 

1.  Theory  and  Method 

The  present  calculations  of  Hugoniot  temperatures  follow  the  previous 
method  with  less  restrictive  assumptions;  only  the  principal  differences 
will  be  discussed  here. 

We  assume  that  the  Hugoniot  equation  of  state  for  C-7  epoxy  given  by 
Eq.  (B-l)  and  that  the  Mie-Gr&neisen  equation  of  state 

V  •  <B-4> 

which  relates  the  Hugoniot  data  (subscript  H )  to  isentrope  data  (subscript 
JO  at  a  given  volume  V,  is  valid  for  this  solid.  The  isentropes  PK  (V)  are 
constructed  by  an  iterative  procedure  from  Hugoniot  data  as  before.  In 

3  J.  M.  tslsh  nd  It  H.  Christian,  Rhys.  R».  97,  1544  (1955). 


75 


this  case  we  have  an  experimental  value  at  one  point  for  T  such  that 
T  =  1.76  * 20 %  at  v  =  0.63.  This  value  was  measured  by  using  the  reflected 
shock  wave  technique*  and  assumes  the  solid  obeys  the  Mie-Griineisen  equation 
of  state.  Also,  since  we  have  a  value  for  T  at  initial  conditions,  deter¬ 
mined  from  elastic  and  thermal  data  for  Lucite,  and  another  determined 
partly  from  C-7  data,  we  assume  straight-line  P  equations  as  functions  of 
relative  volume  between  experimental  and  initial-conditions  points  (Re¬ 
gion  A).  Two  of  these  curves  pass  through  the  upper  and  lower  limits  of 
experimental  probable  error.  Outside  the  range  of  the  experimental  points 
we  assume  T  to  be  constant  (Region  B)  at  one  of  the  three  values,  as  shown 
in  Fig.  B-l.  The  equations  of  the  four  T  curves  and  the  range  of  their 
assumed  validity  are  as  follows:5 


ri  “  P0*U.30  “  3.  30 v )  for  1.00  >  v  >  0.63; 

rt  =  1.76  for  0.63  >  v  >  0.49  . 

°S 

r,  *  r,  («/)  -  r0J(  — —  u  -  v)  for  1.00  >  »  >  o.63; 

2  1  0.292 

r2  -  1.408  for  0.63  >  v  >  0.49  . 

°2 

P.  *  r.  (v)  +  r'  - - (I  -  v)  for  1.00  >  v  >  0.63; 

3  *  0. 38 

Tj  *  1.866  for  0.63>v>0.49 


(B-5) 


(B-6) 


(B-7 ) 


r4  -  (2.94  -  1. 94v)  for  1.00  >  v  >  0.63; 

r4  =  1.76  for  0.63  >  v  >  0.49  . 


In  each  case  the  first  equation  is  for  Region  A  and  the  second  is  for 
Region  6. 

The  results  contain  a  set  of  isentropes,  one  for  each  T  equation, 
and  a  set  of  curves  of  Hugoniot  temperatures  vs.  pressure. 


*  See  text.  Section  IV  B  for  experimental  technique  end  Appendix  C  for  e  celcoletion  of  P. 

5  The  notation  ia  to  be  taken  no  that  P  (»)  ore  an  a  "the  ith  gene  a  an  a  function  of  relative  volune.  ” 


FIG.  B-l  GRUNEISEN’S  RATIO  IN  TWO  REGIONS.  Four  GrJneison’s  ratios 
used  for  regions  A(1.0  ^  v  >  0.63)  and  B(0.63  >  v  >  0.49).  Tho  spread 
in  region  A  at  v  -  1.0  is  duo  to  tho  difforonco  botwaon  initial  data 
of  Lucito  (upper)  and  C-7  (lower).  Tho  spread  in  region  B  is  due  to 
the  upper  and  lower  limits  of  probable  error  of  measurement. 


In  this  calculation,  the  isentrope  temperatures,  TK  ,  are  calculated 
according  to  the  formula 


‘o  a 


=  c  *.  i  =  1,  2.  3,  4  . 


(B-9) 


The  I\<0>  are  the  integrals  containing  and  are  given  by 

.•  r/a,(v)dw 


I‘a)(e)  = 


■  I  “ 


,  i  ■  1,  2,  3,  4;  a  =  A,  B  .  (B-10) 
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Explicitly,  in  the  two  regions  A  and  B,  where  A  corresponds  to 
1  >  v  >  0.63  and  B  corresponds  to  0.63  >  v  >  0.49,  the  integrals  can  be 
shown  to  be6 


I<A)  =  ^(4.30  lnv  -  3.30  +  3.30)  , 

Jj(  B  5  =  IjU)( 0.63)  +  1.76  [lnv  -  n(0.63)] 


(B- 11) 


I^Hv)  ~  (lnv  -  v  +  1)  , 

1  0.292 


I2lB)  =  J2U)(0.63)  +  1. 408 [lnv  -  ln(0.63)]  , 


(B- 12) 


a2 


t(A)  .  jy)(v)  +  — —  ( lnv  "  v  +  1)  , 

3  4  0.38 


I3(b)  *  I3(  4  *(0.63)  +  1.866  [lnv  -  ln(0.63)]  , 


( B- 1 3 ) 


J4U>  *  17,(2.94  lnv  -  1.94v  +  1.94)  , 

ljB)  .  jM)(0.63)  +  1.866[nv  -  ln(0.63)] 


(B- 14) 


In  order  to  obtain  the  Hugoniot  temperatures  Tg  from  the  isentrope 
temperatures  TK,  we  use  specific  heat  data  and  integrate  at  constant  vol¬ 
ume  for  each  value  of  v  from  the  isentrope  to  the  Hugoniot.  That  is, 
from  the  first  law  of  thermodynamics 


dE  =  TdS  ~  PdV 


(B- 15) 


and  the  identity 


TdS  =  CydT  +  T 


(B-16) 


^  Th*  functional  notation  of  the  prtfioni  footnota  holds  for  /.  hara.  Also,  Jj  (0.63)  it ana 
"integral  !{■*)  eralMitd  at  »  =  0.63." 
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we  get  (at  constant  volume) 


dE  =  CvdT  .  (B-17) 

Integrating  from  isentrope  to  Hugoniot  and-  substituting  in  the  internal 
energy  difference  from  Eq.  (B-4)  we  get 

f  Th  ^ on 

I  CydT  =  E„  ~  Ek  =  —  v(P„  -  PK)  .  (B-18) 

\  * 

At  this  point  we  introduce  specific  heat  as  a  function  of  tempera¬ 
ture  as  shown  in  Fig.  B-2.  This  curve  is  obtained  partly  from  experiment 
and  partly  from  theory  as  follows:  We  have  a  first-degree,  least- squares 
fit  of  experimental  data,  for  specific  heat  at  constant  pressure  which  is 
good  in  the  temperature  range  323°K  to  S23°K,  and  is  given  by 

Cp  =  q,'o(<*  +  0T)  ,  (B- 19) 

where  <*=0.331  and  >8  *  1.805-*  10-3/°K.  However,  since  we  are  dealing 
with  temperatures  considerably  higher  than  this  range,  and  since  we  lack 
experimental  data  at  high  temperatures,  we  assume  the  specific  heat  curve 
reaches  the  Dulong-Petit  value  of  6  cal/atom-degree  for  solids  (which  oc¬ 
curs  at  902UK)  and  levels  off  at  this  value  for  the  region  of  interest. 
This,  of  course,  neglects  possible  electronic  contribution  ~at  high~t"e®- 
peratures.  Further,  since  Cy  does  not  differ  greatly  from  Cf,7  we  define 
C^in  two  regions  to  be  (see  Fig.  B-2) 

Region  I:  Cy,  ■  C;„  (<*  +  /8T),  300°K  <  T  <  902°K  ,  (B- 20) 

Region  lit  Cy,,  *  3. 63  *  10' 2  kbar  cm  Vg°K,  902 °K  <  T.  (B-21) 

The  value  in  (B-21)  is  just  the  Dulong-Petit  value  in  the  system  of  units 
we  use  here.  Substituting  these  equations  into  the  left-hand  side  of 
Eq.  (B-18)  we  obtain  equations  for  Hugoniot  temperatures  in  the  two  regions: 

^  J.  C  Slitir,  /n tra Action  to  Oktmieml  Phytict,  McGrnvHill,  lac.,  Nn  York,  1939,  p.  237. 


79 


FIG.  B-2  SPECIFIC  HEAT  IN  TWO  REGIONS.  In  region  1(300  <  T  <  902)  Cy.  it  the 

first-degree  fit  of  experimental  data  oxtrapolatod  to  tho  Dulong-Potit  value  for  C-7. 
In  region  II  (902  <  T)  Cyt,  ii  taken  to  be  constant  at  the  Dulong-Petit  value. 

The  difference  Cp  -  Cy  is  neglected. 


Region  I: 


Tj 


2  V0B 

m 


V 

r. 

i 


(p„  -pk)  +  t«k 


Region  II:  Tg  »  Tg  +  (P„  -  PK) 

Cyil  * 


«  0  . 

(B- 22) 
( B- 23) 


These  are  the  equations  that  are  solved  for  the  Hugoniot  temperatures 
for  each  T.  . 

X 

Some  of  the  computational  results  are  reproduced  in  graphic  form  in 
Fig.  B-3  as  a  family  of  isentropes  (one  for  each  T  used)  along  with  the 
Hugoniot,  Eq.  (B-l),  for  comparison.  The  isentropes  which  develop  a posi 
tive  slope  dPg  / dv  >  0  clearly  become  nonphysical  at  very  high  pressures. 
This  result  is  not  inconsistent  with  the  equation  for  the  slope  of  the 
isentropes,  which  contains  terms  dependent  upon  P  and  dPg/dv  in  such  a 
way  that  a  reduction  in  the  magnitude  of  either  or  both  of  these  quanti¬ 
ties  produces  well-behaved  isentropes.  Furthermore,  the  fact  that  this 
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occurs  where  we  lack  experimental  data  on  T  and  PH  implies  that  the  assump¬ 
tions  on  r  and  the  extrapolation  of  the  curve  beyond  the  experimentally 
determined  region  may  be  invalid  at  very  high  pressures.  Further  investi¬ 
gation  is  needed  in  the  high-pressure  region.  We  note  that  the  isentropes 
are  quite  sensitive  to  variations  in  T  at  high  pressures. 

Further  intermediate  results,  the  temperatures  TK  along  the  isentropes, 
are  given  in  Fig.  B-4.  They  are  labeled  according  to  which  T  was  used  in 
their  calculation. 


FIG.  B-4  AUXILIARY  COMPUTATIONAL  RESULTS.  Temperatures  along  the  various 
isentropes  of  C-7  labeled  by  the  T  used  in  each  calculation. 
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the  principal  results,  the  Hugoniot  temperatures  at  very  high  pres¬ 
sures,  are  given  in  Fig.  B-5.  These  Hugoniot  temperature  curves  do  not 
exhibit  much  sensitivity  to  variations  in  T  as  compared  to  that  of  the 
isentropes  (Fig.  B-3).  From  an  examination  of  the  isentropes  (Fig.  B- 3) 
it  appears  that  the  best  or  most  physically  valid  P  used  is  T2  ,  which  im¬ 
plies  that  the  most  reasonable  Hugoniot- temperature  curve  is  that  labeled 
T2  in  Fig.  B-5.  The  breaks  in  the  curves  of  Fig.  B-5  at  approximately 

O 

900  K  are  a  peculiarity  of  the  computer  program,  which  did  not  smoothly 
join  regions  I  and  II  of  the  specific  heat  data. 

Figure  B-6  compares  the  present  7J,  calculations  with  the  previous 
ones2  made  for  Lucite  up  to  125  kbar.  The  present  temperatures  are  seen 
to  be  somewhat  lower  than  those  calculated  previously;  this  is  in  large 
part  because  of  the  higher  values  of  specific  heat  used  for  the  C-7. 


2#  P.  &  Moftvek,  ibid. 
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FIG.  B-5  HUG0NI0T  TEMPERATURES  OF  C-7  EPOXY  FOR  FOUR  DIFFERENT 
GRUNEISEN'S  RATIOS.  Regions  I  and  II  are  the  temperoture  domains  of  the 
two  different  specific  heat  expressions.  See  Figs.  B*1  and  B-2. 


FIG.  B-6  HUGONIOT  TEMPERATURES,  LOW  PRESSURE  REGION.  Comparison  of  tho 
present  C-7  temperature  calculations  up  to  125  kbar  (solid  curves)  with  those  done 
previously  with  strictly  Lucite  data  (dashed  curves)  assuming  Cy  constant. 
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APPENDIX  C 

CALCULATION  OF  GIHJNEISEN’  S  RATIO,  I\ 
FROM  REFLECTED  WAVE  DATA 


| 

1 

? 

f 
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APPENDIX  C 

CALCULATION  OF  GRDNEISEN’  S  RATIO,  T, 

FROM  REFLECTED  WAVE  DATA 

The  Griineisen  ratio,  T,  expresses  the  relation  of  pressure,  P,  to 
internal  energy,  E,  of  two  states  P^Ey  and  P2E2  at  volume,  V,  of  a 
thermodynamic  system; 1 


Pl  -  P2  -  y  <*1  -  V 


(C-l ) 


Thus,  if  a  material  is  compressed  to  the  volume,  V,  at  two  different 
Hugoniot  states  PyVE j  and  P2VE2,  T  can  be  measured.  The  reflected  wave 
shots  (see  Section  IV  of  test)  accomplish  state  P2VE2  by  compressing  C-7 
from  an  initially-shocked  state  as  shown  in  Fig.  C-l.  The  state  PyVyEy 
is  reached  by  one  shock  from  initial 
state  P0V0£0,  whereas  state 
p2V2E2(Vy  =  V2)  is  reached  from 
shocked  state  PLVLEL;  i.e 


the 


reflected  wave  Hugoniot  is  compared 
to  the  direct  shock  Hugoniot.  The 
value  of  r  is  then  determined  from 
Eq.  (C-l)  as 

(Pl-P,) 


where 


<*i  -e2) 


+  TPi(K«  "  V 


FIG.  C-l  PRESSURE-VOLUME 

DIAGRAM  OF  C-7  STATES, 
REFLECTED  WAVE  SHOTS 


1*  Nw  H*  Ric«»  Ra  6#  IM^«CRt  aid  Ja  Ha  Vilib,  Solid  Sfte  Phytic* ,  Vol  a  6*  Fa  S«itt  ud  Da  Tanibvll, 
tdla,  Actdraic  PrtSSf  Ra  Ta,  1958,  pa  40a 
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and 


£0  +  2  (P2  +  P°)(V0  ~  V2} 

*0  +TP0(V0  -  +T  (P2  +  P0)(Ki  "  V'2) 


Since 


V ,  =  V, 


V, 


r  =  nPj  -  Jyjj V? o  *  K)  "7po(vo  -  "  j  <p2  +  "  K>] 

H^)'p;(i'9-<p’+',;,(f  0]  (C 

(L.- 


2K 
*0 


The  pressure  and  volume  in  state  PqV '0E'a  are  calculated  from  experi¬ 
mentally  determined  free  surface  velocities  in  the  usual  manner.  The 
pressure  in  state  P^ V 2 E-^  is  assumed  equal  to  -that  in  the  reflecting 
plate  which  is  calculated  from  measured  plate  shock  velocity,  free  sur¬ 
face  velocity,  and  initial  density.  The  volume,  V2  ■  l/p2,  is  calculated 
from  the  Rankine-Hugoniot  conservation  equations  which  apply  across  the 
reflected  shock  front: 1 


Pi(t/2  +  «;>  -  p2(u j  +  uj)  (e-si 

and 

P2  -  P'o  -  p'0(U2  +  uj  )(uj  -  U2)  (C-4) 

where  pj,  Uq  are  the  density  and  particle  velocity,  respectively,  in  state 
PjKgfiJ  and  pi ,  U2 ,  and  u2  are  the  density,  shock  velocity,  and  particle 
velocity,  respectively,  in  state  P2^2®2-  The  density  p2  can  be  found  by 
solving  Eq.  (C-4)  for  U2  and  substituting  in  Eq.  (C-3). 

A  value  of  T  can  be  calculated  as  outlined  above  from  shot  No.  9921 
by  first  calculating  the  specific  volume  K2  at  pressure  P2  according  to 

Eq,  (C-3)  and  values  of: 
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P'g  =  59  kbar 

V'0  =  0.63  ccVg 

u'q  =  1.14  mm/fj-sec 

Pj  =  1^5  kbar 

u 2  =  0.36  mm/fisec 


From  Eq.  (C-3): 


U ,  *  3.73  mm//^8ec 


which,  upon  substitution  in  Eq.  (C-4),  gives  a  specific  density  of 

P2  m  77“  =  1-85  g/cc 


Since  Fj  •  Fj ,  a  value  of  Pl  can  be  obtained  from  existing 
Hugoniot  data. 

Substitution  in  Eq.  (C-2)  then  yields  a  value  of  T  «  1.76  at 

F/F#  -  0.63. 
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